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ABSTRACT 

The flat spectrum radio quasar (FSRQ) 4C -1-21.35 (PKS 12224-216) displays prominent nuclear 
infrared emission from « 1200 K dust. A 70 - 400 GeV flare with ~ 10 min variations during half 
an hour of observations was found by the MAGIC telescopes, and GeV variabiUty was observed on 
sub-day timescales with the Large Area Telescope on Fermi. We examine 4C -t-21.35, assuming that it 
is a source of uhra-high energy cosmic rays (UHECRs) . UHECR proton acceleration in the inner jet 
powers a neutral beam of neutrinos, neutrons and 7 rays from p^ photopion production. The radiative 
efficiency and production spectra of neutrals formed through photohadronic processes with isotropic 
external target photons of the broad line region and torus are calculated. Secondary radiations made 
by this process have a beaming factor ex (5^, where fo is the Doppler factor. The pair-production 
optical depth for 7 rays and the photopion efficiency for UHECR neutrons as they pass through 
external isotropic radiation fields are calculated. If target photons come from the broad line region 
and dust torus, large Doppler factors, (5d ^ 100 are required to produce rapidly variable secondary 
radiation with isotropic luminosity > 10^^ erg s~^ at the pc scale. The 7-ray spectra from leptonic 
secondaries are calculated from cascades initiated by the UHECR neutron beam at the pc-scale region 
and fit to the flaring spectrum of 4C -t-21.35. Detection of > 100 TeV neutrinos from 4C -1-21.35 or 
other VHE blazars with IceCube or KM3NeT would confirm this scenario. 
Subject headings: black holes: jets — gamma rays: theory — radiation mechanisms: nonthermal 



1. INTRODUCTION 

Three flat sp e ctrum radio qua sars, 3C 279 
( Albert et al.l[2008l: lAlek sic et al."2011a1. PKS 1510-089 
( Wagner et all 120101: iCortina et al. 2012), and PKS 
1222-H216 (jAleksic et al.ll2011bD . have been detected in 
the very-high energy (VHE; > 100 GeV) regime. This 
last object, a blazar at redshift z = 0.432 also known 
as 4C -1-21.35, is remarkable for day-scale GeV flares 
monitored with the Large Area Telescope (LAT) on 
the Fermi Gamma-ray Space Telescope in 2010 April 
and June (iTanak a et al.l 120111 ) t hat coincided with a 
period of high infrared activity (jCarrasco et al.l I2010D . 
Moreover, the Major Atmospheric Gamma-ray Imaging 
Cherenkov (MAGIC) telescopes observed a VHE flaring 
episode during a period of observation lasting about 
one-half hour on 2010 June 17 that showe d R^ factor-of- 
two va riations on a 10 minute time scale (lAleksic et al.l 
l2011bf) . The rapidly variable 7 rays reached large 
apparent isotropic 7-ray luminosities Lj > lO^'' erg s~^ 
in the VHE band and > 10^^ erg s'^ above 100 MeV. 

The intense and highly variable 7-ray fluxes measured 
from blazars, including 4C -f21.35, demonstrate that 
they are powerful accelerators of energetic particles. En- 
ergetic leptons make intense fluxes of synchrotron emis- 
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sion, with associated 7 rays formed by Compton pro- 
cesses from t he same electrons. As sh own by Fermi LAT 
observations (jAckermann et an2011f ). the two-year time- 
averaged values of L^ of dozens of FSRQs exceed 10^^ erg 
s~^. During flaring episodes, values of Lj > 10*^ erg s""'^ 
are not unusual (the record-setting value is L-y « 2 x 10^° 
erg s~^ from the N ovember 2010 outburst of 3C 454.3; 
lAbdo et al.ll20lial ). In addition to relativistic leptons, 
hadrons are also likely to be accelerated in blazar jets 
and make secondary 7 radiations by interacting with tar- 
get photons of the internal and surrounding radiation 
flelds. This emission will contribute to the formation 
of th e spectral energy distribution (SED) (see iBottchei] 
120101 for a recent review), and the accelerated hadrons 
or secondary neutrons could escape from the acceleration 
region to become ultra-high energy cosmic rays (UHE- 
CRs) with energies > 10^^ eV. 

Typical quasars have broad optical and ultraviolet 
emission lines that originate from the broad line region 
(BLR) near the supermassive black hole. The rapid vari- 
ability, as short as a few hours for 4C +21.35 as found 
in the Fermi-L AT data (|Foschini et al.ll2011a irbl). indi- 
cates that the high-energy 7 rays are formed close to 
the black hole within the BLR. But if this is the case, 
such high-energy photons would be strongly attenuated 
through the 77 — >■ e+e" process against scattered line 
and accretion-disk radiation at the sub-pc scale. We ar- 
gue that the contradictory behaviors of rapid variability, 
large luminosity, and a production region distant from 
the central nucleus can be reconciled if the inner jet accel- 
erates UHECRs that escape as neutrons, which then un- 
dergo photo-hadronic breakup in the infrared (IR) field 
of the dusty torus to make rapidly variable synchrotron 
7 rays at GeV - TeV energies on the pc scale. Further- 



more, UHE 7 rays made in the inner jet can contribute 
to the cascade at the pc scale. 

The recent observations of 4C +21.35 are summarized 
in Section[2] An overview of the model, with estimates of 
the radiative efficiency and secondary photohadronic pro- 
duction spectra, is given in Section [31 A detailed treat- 
ment of secondary production of neutrals is presented 
in Section |4l giving beaming factors and efficiencies of 
secondaries made in relativistic jets. Production spectra 
of neutrons and neutrinos by accelerated UHECRs are 
calculated. In Section \5\ formulas for and calculations 
of pair-production opacity and photohadronic efficiency 
of 7 rays and neutrons as they travel through the BLR 
are given. We show at the end of this Section that the 
BLR is highly opaque to GeV and TeV 7 rays, while 
TeV photons are still strongly attenuated at the pc scale. 
Ultra-high energy protons and neutrons escaping from 
the inner jet are shown to lose significant energy through 
photopion production at the pc scale. In Section [51 the 
parsec-scale magnetic field values are derived that allow 
for sub- 10 minute variability of lepton synchrotron radi- 
ation from electrons and positrons (hereafter referred to 
as "electrons") formed as secondaries of photohadronic 
neutron production. Numerical fits to the MAGIC data 
of 4C -f21.35 are presented. Discussion and summary 
are given in Section [T] 

We make the simplifying assumption throughout that 
the radiation fields are locally isotropic. This assump- 
tion becomes increasingly poor if the radiation energy 
density is dominated by anisotropically distributed pho- 
tons impinging from behind. Thus, provided the emit- 
ting region is within the BLR radius (i?BLR) for processes 
involving atomic-line or reprocessed accretion-disk radi- 
ation, or within the pc scale for processes involving the 
IR radiation field of the dusty torus, these results give an 
accurate determination of the 77 opacity and photopion 
efficiency. Note also, that for 4C -^21.35 at z = 0.432, 
the luminosity distance dL = 2335 Mpc = 7.20 x lO^^ 
cm in a ACDM cosmology, where the Hubble constant 
h = 0.72, in units of 100 km s^^ Mpc"\ f7,„ = 0.27 is 
the ratio of matter density to the critical density, and 
r^A = 0.73 = 1 — ilm is the ratio of cosmological constant 
to the critical density. 

2. OBSERVATIONS OF 4C -1-21.35 

In this section, the relevant observations and radiation 
fields in the nuclear region of 4C -1-21.35 are described. 

2.1. Observations 

The discovery of VHE emission betw een ^ 70 and 
400 GeV wit h the MAGI C telescopes (jMariotti et al.l 
I2010t JAlcksic et al.l I2011bf ) was anticipated by the re- 
port of VHE photons found in Fermi LAT data 
during a high state of 4C -1- 21.35 in April 2010 
(|Neronov. Semikoz fc Vovk1l2010[ ). The source was tar- 
geted with the MAGIC telescopes from 3 May 2010. 
The spectrum of the major VHE flare observed with 
MAGIC on 17 June 2010, which observed for only - 30 
minutes, is well described by a single power law with 
photon number index ar — 3.75 that is hardened by 
about one unit whe n corrected for attenu ation with a 
low-intensity model (jKneiske fc Dolel[2010D of the extra- 
galactic background light (EBL) at infrared and optical 
frequencies. The remarkably short 10 minute variability 



(in comparison, 3C 279 displ ayed weak VHE v ariabil- 
ity on time scales of days; see I Albert et al.|[2008[ ) points 
to a production site near the central supermassive black 
hole where coUimation by extreme processes can occur 
(jTavecchio et al.ll2010D . Variability associated with the 
dynamical timescale of a lO^Mg Solar mass black hole 
is idyn = (1 + z)rsh/c = (1 + z)10''Mg s, where Tgh is 
the Schwarzschild radius of the black hole, whereas vari- 
ability on time scales shorter than tdyn by factors of 10 
or more has been o bserved in the rapidly variable BL 
Lac objects Mrk 501 (lAharoman et al.lboO^ PKS 2155- 



304 |Aharoman et al.ll2007[) . and Mrk 421 (jFossati et all 
I2008lf) . For 4C -1-21 .35, with a black h ole mass estimated 
at « 1.5 X 10^ Mq (|Wang et al.ll2004D . the corresponding 
timescale is « 2000 s, so the VHE measurements indicate 
a marginally hyper- variable source. 

A higher black-hole mass is obtai ned from recent op - 
tical spectroscopy of J1224-t-2122 bv lShaw et all (|2012D . 
using scaling relations derived from reverberation map- 
ping of radio-quiet AGNs. From analysis of the FWHM 
of the H/? line, they deduce a black-hole mass for 4C 
+21.35 of 7.8(±2.3) x IO^Mq, which does not include 
uncertainties in applying scaling relations to blazars. 
Th is value is co n sidera bly higher than the mass reported 
by iWang et"al] (|2004D . and would relax the energetics 
based on the Eddington luminosity, but make the short 
timescale variability more problematic. Here we use the 
lower mass, and return to this point in Section 7. 

The BLR radius -Rblr can be written in 
terms of the accretion-d isk luminosity Ldisk as 
(iGhiseUini fc Tavecchiol [20081 ) 



R 



BLR 



?^10^ 



-t^disk 



2045 erg g-i 



(1) 



Based on t he H/3 luminosity, Lh^ = 2 x 10"''^ erg s 



Tanaka et al. (2011) de rive from the sca ling relations of 
Wang et al.l (|2004i ) and iFan et al.l (poM ) that the BLR 
luminosity I/blr « 5 x 10''^ erg s~-^, and that the ac- 
cretion disk luminosity Ldisk is an order of magnitude 
greater, Ldjsk ,^ 5 x lO'^^ erg s~^. In comparison, 
JTavecchio et al.l (|2011h estimate that Ldisk = 5 x 10"'^ erg 
s~^, assuming the Swift UVOT spectrum from 4C +21.35 
is the optically thick accretion-disk radiation field, in 
which case i?BLR ~ 7x 10^^ cm « 0.2 pc. This represents 
super-Eddington luminosities, given that ^Edd ~ 2 x 10"*^ 
erg s^^ unless the black- hole mass is > 4 x 10^ Mq. We 
use the smaller disk luminosity here, noting that a larger 
Ldisk and BLR radius will improve photohadronic effi- 
ciency. 

More precise relations depend on specific lines and as- 
sumptions about the geometry of th e BLR (see, e.g., the 
recent discussion bv lFoschm]||2011f ). In particular, we 
could use the rest-frame 5100 A or 1350 A continuum, 
or the high-ionization Lyman a line , as seen strongly 
in, e.g., 3C 454.3 (Bonnoh et al.ll2011[ l. For definiteness, 
we focus on the Ly a line to illustrate BLR line effects, 
which typically car ries > 20% of the BLR luminosity 
([Francis et al.|[l99l[ ). Approximating the Ly a/BLR line 
luminosity Li^ya = O.lLdisk- and defining 



-^Lya/Ldisk 
(-RLyaZ-RBLR,)^ 



0.1(j). 



(2) 



gives a Ly a photon field with energy density ULya ~ 
LLya/47ri?2y„c ^ 0.026<^_i erg cm-^. 

The magnetic field B' in the comoving jet frame can 
be determined if the GeV 7 rays are Compton-scattered 
external radiation and the lower energy optical radiation 
is synchrotron radiation made by the same nonthermal 
electrons making the 7 rays at GeV energies. The Comp- 
ton dominance parameter Aq is defined as the ratio of 
the 7-ray luminosity (assumed to originate from exter- 
nal Compton scattering) to the radio/X-ray synchrotron 
luminosity. For 4C +21.35, Aq reaches values as large 
as K, 100 (see multiwavelength SEP using contempora- 



neous data in lTavecchio et aT]l2011l ). If the external ra- 
diation field scattered by these electrons are the Ly a 
photons, then Aq « M Lyg/^B' with v!^^^ « 4r2uLyQ/3, 
and m'b = B'VStt Ce.g. JSikora et al.ll2009H . Thus the co- 
moving magnetic field for a blazar jet with bulk Lorentz 
factor r is 




(3) 



Modeling 5 - 35 /i Spitzer SPSS, 2MASS and S wift 
UVOT data of 4C -^21.35, IMalmrose et"aI1 (I20T1I) de- 
compose its spectrum into a nonthermal power-law and 
two-temperature dust model. Hot dust with T « 1200 
K radiates w 8 x 10"*^ erg s~^ from a pc-scale region, 
and a second warm dust component radiates sa 10^^ erg 
s~^ at T K, 660 K on the same scale. The IR com- 
ponent, which has been proposed as a target in external 
Comp ton scenarios (Blazciowski et al. 2000; S ikora et al.l 
|2009( ). is especially important for UHECR scenarios, not- 
ing that the relative energy densities of the BLR, the 
dust fields, and the CMBR (at z = 0) are, in units of erg 
cm-3, UBLR ^ O.O260_i, Udust = 3 X 10-3(Ljfl/1046 erg 
s^^)/(i?/pc) , and ucmbr = 4 x 10~^^, respectively. We 
use th e same parameters as derived by IMalmrose et al.l 
(|2011f) in our model. 

2.2. Target radiation fields 

Candidate target radiation fields are the cosmic mi- 
crowave background radiation (CMBR) , with dimension- 
less temperature 8cmbr = kBTcMBR./'meC^ — 4.6 x 
10-i°(l + z), 1200 K and 660 K dust (Gdusti = 2 x 
10~^ and 0dust2 — 1-1 X 10"^, respectively), quasi- 
thermal scattered 10 eV disk emission (8acr disk = 2 x 
10~^). For Ly a line radiation — neglecting broadening, 
which is unimportant for these calculations — CLya = 
10.2 eV/m^c'^ = 2 x 10"^ Using the GZK energy 
of w 6 X 10^^ eV as a yardstick where A isobar exci- 
tation by 2.725 K CMBR photons is a source of opac- 
ity and energy loss, the IR dust and accretion disk/Lya 
radiation are effective targets for photopion production 
from cosmic-ray proton and neutron (nucleon) interac- 
tions with Ep > 10^'^ eV and Ep > 10^^ eV, respec- 
tively. The energy densities of the target radiation fields 
are defined by the usual relation uq w L/A-kE?c where R 
is the distance from the black hole0 

* This may underestimate the mean energy density inside a uni- 
form spherical distribution o f emitters, which is a factor 2.24 larger 
HAtovan &: Aharomanll996l ). We use, however, the simpler relation 
throughout. 



The luminosity of photopion secondaries is propor- 
tional to the effective photon number density for pho- 
topion production. For monochromatic line sources, 
rieff = Uq / Ci^meC^ , and for thermal blackbody or gray- 
body radiation fields, 'T.e// = uo/QrrieC'^, as shown be- 
low. Here, e^, is the mean energy of photons in the unit 



of 



rUeC . The effective photon number density for the 



CMBR is 



CMBR 



WCMBR 



"leC^OcMBR 



= 1.1 X lO'Xl + zY cm"-^ (4) 



assuming the temperature of the CMBR at the present 
epoch is 2.725 K. The effective photon density for p7 
scattering of the Ly a field is 



Lycx 



Uhya 



eLyc 



1.6 X 10^ 



(5) 



The BLR radiation field due to the accretion-disk radi- 
ation scattered by BLR gas is, when approximated as a 
graybody, given by 



"BLR 



BLR _ 



acr disk 



1.7 X 10^r_i 

(Oacr disk/2 X 10-5) 



(6) 

assuming an effective Thomson scattering depth of 
0.1t_i. For the two dust radiation fields, the effective 
photon densities for photopion p7 interactions are 



and 



dustl 
^eff 



dust2 



Udustl 



meC^edustl 



Mdust2 
meC^Odusfl 
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3.1 X 10-*L45 _s 
cm 



-'•'pc 



(7) 



(8) 



Here the luminosity of the IR radiation, assumed to be 
radiated on a size scale of Rpc pc, is denoted 10-' L^ 
erg s~^. To order of magnitude, rieff = uo/meC^e ~ 
uo/meC^Q ~ 10^ - 10^° cm-^, and the CMBR is negli- 
gible. 

For sufficiently energetic cosmic-ray nucleons, the 
warm dust radiation field provides the densest and most 
important target photon field for photopion production. 
For dissipation of energy on the pc scale, only the highest 
energy nucleons with E > 10^^ eV lose energy effectively 
by photohadronic processes with IR dust photons, be- 
cause the BLR radiation does not extend to the pc scale, 
and there is no other sufficiently strong radiation field at 
the pc scale to extract energy efficiently from these lower 
energy protons. 

3. MODEL AND ESTIMATES 

The model is outlined, followed by a simple estimate 
of photohadronic efficiency and a simple derivation of 
secondary neutral production spectra made by photo- 
hadronic interactions of cosmic-ray jet protons in the in- 
ner jet. 

3.1. Model description 

This model is motivated by the hypothesis that that 
bla zars and radio galaxies accelerate UHECRs (e.g. , 
;Mannheim &: Biermann. ,1992; .Berezinskv et al.l 120061 : 



-f UHECR 




IR Dust Torus 



Accretion Disk 

Fig. 1. — Cartoon of the system. Ultra-high energy neutrals, 
made through photohadronic production of cosmic-ray protons ac- 
celerated in the inner jet of 4C -f21.35, escape to form a beam of 
neutrons, neutrinos, and 7 rays. Subsequent photohadronic and 77 
interaction of the neutrons and UHE 7 rays with IR photons from 
the dust torus induce a high-energy lepton beam that quickly radi- 
ates secondary VHE photons to make the variable 7-ray emission 
in 4C -1-21.35. UHE neutrons and ions escaping from the system 
become UHECRs, and the jet radiation is accompanied by an as- 
sociated flux of UHE neutrinos. 



iDermer et all 120091: iMurase fc Takamil I2OO90 ; confirma- 
tion of this hypothesis is the detection of UHE neutri- 
nos during blazar flares. Electromagnetic signatures can 
also provide crucial evidence in support of UHECR pro- 
duction in blazars, and the presence of distinct hadronic 
emission signatures in the variable 7-ray SEDs is con- 
sistent w it h the SED sha pe but not u nambiguous (see 
lB5ttcheH (|2010l) or, e.g., lAbdo et all [20 11b , for Mrk 
421). Detection of a multi-TeV radiation signature 
in extreme high-synchrotron-p eaked blazars like lES 
0229+200 (jMurase et all l2012| ) at energies where the 
EBL should suppress this emission would support emis- 
sion generated by UHECR protons accelerated by blazars 
(|Essev fc Kusenkoll2010HEssev et al.ll2010Ll201lD . 

Fig.[l]shows a cartoon illustrating the underlying basis 
of this model. UHECR protons and ions are accelerated 
in the outflowing relativistic plasma formed in the inner 
jet. Interactions with internal synchrotron and exter- 
nal quasi-isotropic radiation field cause ultra-relativistic 
protons accelerated in the outflowing jet to undergo pho- 
topion losses, with the production of escaping n eutrons, 
7 rays, and neutrinos (jAtovan fc Dermerl 120031 ) . Subse- 
quent interactions of the escaping neutrals with IR torus 
photons on the pc scale make highly beamed leptons pro- 
duced as secondaries of 77 interactions of UHE photons, 
and as secondaries of photopion interactions of UHE neu- 
trons. The relativistic leptons make beamed synchrotron 
and Compton scattered radiation. These emissions, be- 
ing highly beamed, will vary on the production timescale 
determined by the behavior in the inner jet. Off-axis 
emission from misdirected leptons will not smear and 
enhance the timescale because it is produced by ultra- 
relativistic electrons that are not significantly deflected, 
so the contribution to the observer originates only from 
the particles traveling directly toward the observer. 

3.2. Radiation efficiency estimates 

We make simple estimates of the radiative efficiency for 
the production of photohadronic secondaries when UHE 
protons interact with photons from isotropic radiation 



fields external to the jet. In our estimates, the energy 
density of the external radiation field is denoted uq and 
the characteristic spatial size scale of the external radia- 
tion field is R. For thermal radiation fields, e^ = 2.708. 

Suppose that a relativistic jet moving with bulk 
Lorentz factor F is filled with cosmic-ray protons and 
ions. The comoving dynamical time for traveling through 
an external radiation held of is t'^ = R/Tc. In the jet 
frame, the energy density of external radiation field is 
Uq = F^mq, and the mean photon energy is e^ = Fe^. 
The threshold Lorentz factor 7' ^-^^J. for photomeson pro- 
duction by protons with comoving Lorentz factor 7' 
interacting with photons with energy e^ is given by 
7'p,thr ~ ethr/eJt, where ethr « m^/me and m^ = 140 
MeV is the pion mass. 

The energy-loss rate of ultra-relativistic protons 
through photopion processes is therefore 

/'~W ^ V ^ ~ (-^P70"P7)'^"0 j^ ,'7,- \ p/ 

>7 \lp > 7 p.thrJ ~ 1 I ~ (,-".p7Crp7JCi (^ 



Uq 



TleC^e^' 



(9) 

where a = Kp^Up-j = 70 /ib is the product of the in- 
e lasticity and the cro s s section for photomeso n losses 
(Beg elman et al.l 119901 : lAtovan fc Permed I2003D . Note 
that the energy-loss timescale is proportional to the ef- 
fective photon densi ty Ueff = uo/meC ^e-i,. Employing 
the approximation of lAtovan fc Dermeii ([2003) here and 
in the following, we take ethr = 390. The photopion pro- 
duction efficiency is therefore simply 



Vp-/il'p) = t'dyjt'p^ilp) ~ O-i? 



uo 



meC^ei, 



(10) 



provided 7' > 7' j-j^j.. The Lorentz factor of protons 
measured in the stationary frame of the black-hole jet 
system (as if they had escaped the jet) is 7 = F7', so the 
threshold Lorentz factor is 7p,thr — ^thr/e*- 

Following Equations ([5]) - ([8|), we write ?^e//,g — 
(uo/r7ieC^e*)q = [{rieff/lO'^ cm^'^)]. For interactions 
in the BLR with Ly a photons or scattered accretion 
disk radiation, r]p^{'yp) sa 0.007ne//,9i?i7 for escaping 
protons with energy Ep > 2 x 10^^ eV. Here, R = 
10^^i?i7 cm. For interactions with IR photons from the 
dust torus, r]pj{'yp) w 0.07ne//,io^i7 for protons with 
^p ^ 7 X 10"'^'' eV. For larger dissipation radii, one has 
Vpiilp) ~ '2neff, loRpc, which means that UHECR pro- 
tons in the jet plasma can dissipate essentially all of their 
energy into secondaries while traveling through the radi- 
ation field of the dust torus. But note that in this case, 
the UHE proton beam has to be maintained throughout 
the entire pc-scale length of the external radiation fleld, 
in which case the emission would vary on timescales no 
shorter than - R/T'^c ~ lO'* si?pc(F/100)"^ 

Thus we see that the efflciency to produce photomeson 
secondaries through interactions of UHECR jet protons 
with photons of an external radiation field is several % in 
the inner jet, counting both line and scattered ph otons, 
and noting that Ldisk ~ 5 x lO'*^ erg s~^ (Tanaka et al.l 
l2011f ). so that i?BLR ~ 2 X 10^^ cm, from Equation 
([!]). This holds for > 10^'^ eV protons. Under opti- 
mistic conditions, the photopion efficiency of cosmic-ray 
protons in a blazar jet is of order unity in interactions 
with the IR dust photons, but this applies only to pro- 



tons with escaping energies > 10^^ eV. Note furthermore 
that this is a conservative estimate of photohadronic ef- 
ficiency, because it does not take into account internal 
synchrotron photons, or photons from a sheath-spine 
structure or from d ifferent porti ons of the jet (compare 
IGhiselhni et all [2005 : Georganopoulos fc Kazanaj 120031 
for TeV blazars, though similarly structured jets could 
be found in FSRQs) that could significantly enhance the 
target radiation field and the photohadronic efficiency. 

It is interesting that all the F factor dependencies in 
Equation ()10|) drop out. This shows that we can equiv- 
alcntly calculate the photohadronic efficiency for a pro- 
ton traversing the volume of the external radiation field 
without transforming to and from the comoving frame. 
We can furthermore make a simple derivation of the sec- 
ondary neutron production by only considering interac- 
tions in the frame of the black hole and external radiation 
fields. 

With a photopion energy-loss cross section of ~ 70 /ib, 
the column density of photons above threshold required 
for unity radiation efficiency is A/" « l/{Kp^ap^) ^ 1.4 x 
10^^ cm~^. The column density of Lya photons in the 

BLR is NbLR ^ rieffXyaR ^ W^^Ueff^gRn cui"^ . The 

column density of IR photons from the pc-scale warm 
dust torus is Afdusti ~ rieff^dustiR ~ 3 x lO'^^Ueff^wRpc 
cm^^. The IR photon field of the warm dust clearly 
presents the largest column density for photohadronic 
production, for the most energetic nucleons above the 
photopion reaction threshold. 



3.3. Simple Derivation of Secondary Neutron 
Production Spectra 

The estimate above implicitly assumes that when pro- 
tons undergo photopion interactions, they remain pro- 
tons. In about 50% of the time, however, neutrons are 
formed which escape from the jetted plasma in a coUi- 
mated outfiow. The escaping neutrons deposit energy 
throughout the length of the jet by subsequent photo- 
hadronic interactions. We now estimate the energy loss 
rate of protons into neutrons via photopion production, 
and use that result to make a simple derivation of the 
secondary neutron production spectrum, anticipating the 
more detailed derivation in the next section. 

The inclusive energy-loss rate of protons into neutrons 
only is given by 



dE, 



dt ^P''- 



pxCcJpiEp , 



(11) 



where p is the cross-section ratio, x is the fractional en- 
ergy of the incident proton that is deposited into sec- 
ondary neutrons, and C the neutron multiplicity. For 
neutron production near threshold, p = 0.65, x — 0-8, 
C = 1/2, and the maximum photopion cross section is 
CTp^ — 520±30/ib. Here the proton energy Ep = nipC'^jp, 
and rup = to„. 

The rate for a primary proton to lose energy into sec- 
ondary neutrons is therefore simply 



.(7p) 



I dEp I 

En 



denph(e) / dp{l - p)p{er)x{^T)C,{er)ap^{er) 



c 
2 

^ (12) 

where e^ = 7pe(l — m) is the invariant interaction energy. 
Solving for a monochromatic radiation source 7iph(e) — 
UoS{e - e*)/(meC^e) gives 

dEp^ 



C i Cp^ ^p 



m 



dt 

I ethr 



V^Tpf* 



H 



Ip 



ethr 

2^ 



(13) 
where T = pxCO 

The inclusive energy loss by primary protons with 
Lorentz factor 7^ that is transformed into secondary neu- 
trons with Lorentz factor 7„ is conserved; therefore 



dip 



d£r> 



\ dtdjp 



P'y—^71 



din 



d£n 

dtdin 



(14) 



and ip = Inix- The neutron production luminosity is 



lp{d£p/dtdip)p 



For a stationary-frame 



distribution in the effective volume (where the retarded 
time is considered), Np{ip,^p) of protons differential in 
Ip and direction fip, the secondary production spectrum 
of neutrons, assuming that the neutrons have the same 
direction as the original protons, is 



lnLn{ln:^n) 



— cTa. 



Pi—lpNpiip^ilp) 



ethr 

27pe* 



H 



Ip 



Cthr 

27pe* 



(15) 



Us ing the method of IGeorganopoulos et al.l (|2001l 
120041 ) ■ we transform the comoving jet frame proton spec- 
trum iVp(7p, Q.'p) to the stationary frame jet proton spec- 
trum Npdp, flp) through the relation 



NpiipMp) = 



_ dj\fp{ip,np) _ 3 



dipdQp 



^sI,n;,{i;,%) 



iWp) 



(16) 

where the last relation is a consequence of the assumption 
of proton isotropy in the proper fiuid frame, noting that 

Ip = In/x = Sul'p- Thus 

rUp capjUo 



47r7„L„(7„,Sl„) 



27ne* 



e* 



H 



'rsl,[i'p'N;i%)] 



2e. 



(17) 



Note the (Sq beaming factor arising from the transfor- 
mation of the proton spectrum, equation ()16p . and one 
power each from energy and time. 

From Equation (J17p . we can determine the absolute 
amount of energy in protons, £p = T£' that, while travel- 
ing through a radiation field characterized by UQ/e-^meC^, 

^ The Heaviside function with a single argument is defined by 
H{x) = 1 if a: > and H{x) = otherwise. The second version of 
the Heaviside function is a step function defined by H{x; a,b) = 1 
if a < X < b, and H{x; a,b) = otherwise. 



6 



will produce an apparent isotropic luminosity in neutrons 
Ln = lO"*® erg s~^. If the spectrum of protons is —2, that 
is having equal energy per decade in protons, and extends 
well above the threshold 7p.thr — ^thi/'^e* for neutron 
production, then the fractional energy is proportional to 
the ratio of two logarithmic factors, which has a value of 
order « 1/2, implying 



1 rup cap^uo^ ^ 5^ _£p_^ ^ ^^48^ 



2 TTlp 



r TOdC^ 



(18) 



There is no fundamental reason why the engine should 
be Eddington-limited during outbursts, but this provides 
a fiducial amount of energy that could be generated by 
the central engine. Therefore we assume that the abso- 
lute nonthermal proton energy f p w 2 x 10'*^ ivar/(l + 
z) = 8.4 X W^ erg, for the 4C +21.35 fiare (tvar ~ 600 
s is the variability timescale of the flare and z — 0.432), 
which implies that 






> 



1.2 X 1G« 



(T/0.1)K/e*)3 



(19) 



so (5d ~ r > 10^ in order to make an apparent luminosity 
> 10''^ erg s^^ in neutrons. 

The r factor limit can be relaxed to 5u > 50 if we 
happen to be viewing exactly down the jet (Su = 2T) or 
if the energy content in protons is not limited by the Ed- 
dington luminosity during this outburst, but is instead 
one or two orders of magnitude larger. Alternately, small 
r is possible if the photomeson production occurs mainly 
by target photons produced in inner jets (see below) . In 
principle, there is no real reason that such large values 
of the Doppler factor or F factor are not allowed. Argu- 
ments based on 77 opacity appl ied to the July /August 
2006 flares from PKS 2155-304 (lAharonian et al.l l2007D 
require Sb > 60 (jBegelman et al.ll2008() . while a full syn- 
chrotron/SSC model flt, including 77 and EBL effects, 
implies 5b > 100 for PKS 2155-304 and S^ > 80 for Mrk 
501 w ith a variability timescale of w 10^ s (' Finke et al.l 
I2OO8D . One-zo ne leptonic models f or the VHE emission 
from 3C 279 (|Albert et al.l [20081: lAleksic et all IMIIal) 
furth ermore require large , F > 100, bulk Lorentz fac- 
tors ()Bottcher et al.|[2009l) . Standard arguments relating 
variability time and location of emission site through the 
expression ivar ^ (1 + z)R/r'^c imply F > 500 for emis- 
sion produced at the pc scale, which already undermines 
naive expectations about F and the size of the emit- 
ting region. Measurements of superluminal motion never 
reach such values, but this is based on radio observa- 
tions, which may be tracking the outflow speed of a radio- 
emitting sheath rather than a more rapidly moving spine. 
The necessity of a radical departure from standard mod- 
els to expla in the 4C -1-21.35 resu l t seems unavoidable in 
ah models (iTavecchio et al.ll2011l:rNalewaiko et~alll2012l: 



iTavecchio et al.ll2012l) : here we construct a model within 
a framework that is testable by neutrino telescopes. 

4. SPECTRAL POWERS AND BEAMING FACTORS FOR 
PHOTOHADRONIC SECONDARIES 

The system is now treated more carefully. A stan- 
dard blazar jet scenario is considered (jAtovan fc Dermeil 
120011) . In the comoving frame defined by the existence of 
a quasi-isotropic particle and randomized magnetic-field 
distribution moving with Lorentz factor F at some angle 



9 with respect to the observer, secondaries are formed as 
a result of photohadronic processes, e.g., neutrons (n), 
neutrinos {v), electrons and positrons (e^), photons (7), 
protons (p), and /3-decay electrons (c^) and neutrinos 
(vp). The charged secondaries may be trapped in the jet 
by its magnetic field, but the neutrals (n, z/, and 7 rays) 
escape, forming a neutral beam with a beaming factor re- 
flecting the Doppler factor of the jet and the production 
kinematics ()Atovan fc Dermeill2003[ ). 

We start with Equation (2.44) of IDermer fc Menoiil 
([2009), but instead of treating electron-photon scatter- 
ing, here we treat secondaries of protons with Lorentz 
factor 7p made in photopion interactions with photons 
with dimensionless energy e = hv/rrieC^. The technique 
is, however, general and can be adapted to other sec- 
ondaries, noting that the dimensionless secondary en- 
ergy e^ in TOeC^ units will be written as a fraction of 
X of the primary proton energy nipC^jp. Thus for sec- 
ondary photons, leptons or neutrinos, Sg = X"^p7p/"^e- 
For secondary neutrons, we can directly write, because 
TO„ « Trip, the neutron production spectrum in terms 
of the secondary neutron Lorentz factor 7„ = X7p. The 
direction-dependent emissivity for the production of neu- 
trons through the process p -I- 7 -^ n -I- A" becomes 



,3„,2 



lnLn{"fn,^n) = mpC 7, 



dn 



X (1 - cost/') / djp Npi'^p, rip) 



de 



u{e,n) 



dflr. 



dap^{e,n,jp,il,p) 
d^nd^n 



(20) 

specialized to ultra-relativistic protons, 7p 3> 1. Here 
ij} is the angle between the directions of the interacting 
proton and photon. The approximation we use for the 
differential cross section for p7 photo-pion processes is a 
sum of step functions, given by 

dcrp-y(£, r?,7p,17p) _ 

d^ndVln 



N 

^P7 X! ^iPi^i^^'^ ^i' ej+l)<^(7n - Xtlp)S{^n - ^p) • (21) 
1=1 

Here the relativistic invariant is e^ — 7pe(l — cos'tp), 
cos?/; — fJ,fip + \/l — M^\/l ~ MpCOs(0— 0p), and arccos/x 

and (j) are the cosine angle and azimuthal angle, respec- 
tively, of the photon or (with subscript "p" ) proton. The 
term pi is the fractional cross section written as a ra- 
tio of the value of the step-function cross section to the 
maximum pj cross section Cp^ = 520 ± 30 iA>. The effec- 
tive multiplicity for segment i between invariant energies 
Ci ^ Er < Q+i is denoted Q, and Xi is the mean energy 
fraction of a secondary neutron compared to the original 
proton energy for interactions associated with segment i. 
Substituting Equation (PT|) into (^0]) and solving gives 



InLnijn.^n) 



m. 



^ /- 

2 Y^ Si Pi 



i=l 



Xz 



dil (1 — cosV") 
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Jo 



de 



u{e,n) 



A^p(i^,fi„)if(e;;Q,e,+i), (22) 

Xi 



where 



ej, = — e(l — cosip) 



and cos ip = fifj,n + \/l — f^^y/^ — fJ-n cos (p, taking (^„ = 
without loss of generahty. 
Using Equation ([TS]). we have 



Tn -^n w n ? ^ ^n J 



JV 



TTlp / 47r 



■s'uiii: 



dpi 



xjdn (l_cos7^)y°° A^^^ iV;(7;jif(e;;e-„e-,+i) , 

(23) 
where 7'p, = 7„/xi(5D. 

We now speciahze to the case of a surrounding ex- 
ternal radiation field that is isotropic in the stationary 
frame of the black- hole/accretion-disk system. Therefore 
u[e, O) = u{e)/A'K, and we can let /x„ = 1 without loss 
of generality. After some straightforward manipulations, 
we obtain 



N 



7n^«(7«,^n) 



e 



Vi+l 

where Tj = PidXi, and 
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uie) 
e3 
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(24) 



Vii+i) 



_ Xi<^i(+i) 

27n 



Or, we have 



'friLnilnMn) 



Vi+l 
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J p^ 1 1 Lp 



^T,[7p2,iVp(7p„r!„)]; 



Vi ''Vi -^ai + l ^ 

(25) 

4.1. External isotropic monochromatic radiation field 

For a monochromatic external radiation field, we write 
u(e) — uoS{e — e*). In this case, 

7„Lf'=(7„,^«)= (^) "^Sl J2^.[%'.Ki%^)]^ 
^ '^ i— 1 



y2 2 



This can also be written as 



(26) 



AT 



4.7ni^ij-(7„,f7„) ^ 2Zn}^^iihi SI, ^T. [7;;^A^;(7;)] 



TOeE* 
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[Hil„ 
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|iJ(7» - If) 



(27) 




Fig. 2. — Functions Io{x) (red; Equation (|33j) and l2{x) (blue; 
Equation l|34|l ) with asymptotes (daslied and dotted lines). Inset 
shows functions on a log scale. 



The single step-function approximation with £2 
takes the form 

47r7„i^„ (7„, i2„) w 1 1 dj; 



mee* 



X [7pi^;(7;i)] 



if? 



H{ln - 



2e* ' 



(28) 



with Ti = XiCiPi, 7pi = 7n/xi'^D, and yi = xi^i/'^ln- 
This confirms the approximate derivation. Equation (ITTI , 
and also shows that secondary production of photo- 
hadronic secondaries has a beaming factor ex (^p. 

4.2. External graybody radiation field 

The graybody radiation field, defined as a blackbody 
radiation field times the graybody factor g giving the 
ratio of energy densities of the radiation field under con- 
sideration to that of a blackbody radiation field with the 
same effective temperature T = rUeC^Q/kB , is expressed 
as a spectral energy density in the form 



Wgb(e;6) = 9- 



\l exp(e/e) 



1 



(29) 



The electron Compton wavelength Ac — h/nieC — 2.42 x 
10^ ^'^ cm. The energy density of blackbody radiation is 



Ubb(0) 



15A^ 



e^ 



(30) 



and g = uo/'"bb, by definition. 
Substituting Equation ([^ into ([M|) and solving gives 
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s!,Y.^^Wp'^N;,{JU)]^ 
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[h{xi) ~ hixt+i) + xfj^_-i^Io{xi+i) - xfloixi)] , (31) 
where 

2;»(+l) = n. Q ' (32) 
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and 



loix) = 
ln[l — cxp(— x)] — > 

hix) ^ 



dw 



L exp(w) - 1 
r -Inx , < a; < 1 

[ exp(— x) , a; >• 1 



(33) 



^ exp(w) - 1 
r 2C(3) - a;V2 , a; < 1 

I (2 + 2a; + a;^) Gxp(-x) a; > 1 



(34) 



where C(3) = 1.202 ... is Riemann's zeta function ({n) 
for n = 3. Functions /o(a;) and l2{x) along with asymp- 
totes are plotted in Figure [5] 
In the one-segment approximation with 62 — > oo, 



47r7„Lf(7„,r!„) 



IScCTp^TTlp Uq 



4 ^Sl,T,[j;iN^ij;,)]Iix), (35) 

where x = Xiei/27n6, and 

r2C(3)+a;2(lna;-l/2) , x<l 

/(X) = l2{x)-X^Io{x) ^ { 

[ (2 + 2.x) cxp(-a;) x > 1 

(36) 
Integrals Ii{x) = J dww^ /[exp{w) — 1] also arise in the 
study of anisotr opic Compto n scattering of stellar black- 
body radiation (jDermer fc B ottchcr 2006). The large- a: 
asymptote is accurate to within 20% even for small x 
(Fig.©. 

4.3. Secondary photopion production cross section 
An N — 2 segmented c ross section was u sed in 



lAtovan fc Dermeii (|2003[) and iDermer fc MenonI (|2009[ ) 
which accounted for A(1232) resonance and multi- 
particle production in p+"f — > 7r + A" interactions. We ex- 
tend the N — 2 segment approximation to iV = 3, which 
now better accounts for the production of heavy reso- 
nances. Parameters are given in Table [T] The number of 
segments in the approximation can be made arbitrarily 
large, but energy dispersion, amounting to another inte- 
gration, is required to reproduce the accuracy of Monte 
Carlo simulations. The minor loss in accuracy is compen- 
sated by code flexibility. Detailed numerical calculations 
of neutrinos in relativistic jet sources, where the multi- 
pion production is taken into account, are found in e.g., 
iMiicke fc Protherod (|200l : iMurasd (|2007l ). 

4.4. Jet power and energetics 

The apparent isotropic jet power during a short flar- 
ing episode of a black-hole jet can be large, and must 
be > lO''* erg s~^ in order to reproduce the value of 
L^ « lO''^ erg s~^ observed in VHE 7-rays from 4C 
+21.35 during the flare (the GeV 7-rays can, and are 
likely to be, produced in the inner jet; see Section 7). For 
a continuous jet, which makes persistent emission over 
long periods of time, the jet power is likely to be smaller, 
but still has to exceed « 2.5 x 10"*^ erg s~^, which is the 
two-year time-averaged 7-ray luminosity measured with 
the Fermi-LAT for 4C -^21.35 (,Ackcrmann ct al.,,2011.) . 



TABLE 1 

N = 3 SEGMENT APPROXIMATION PARAMETERS FOR NEUTRONS AND 

NEUTRINOS FORMED AS SECONDARIES IN PHOTOPION INTERACTIONS 

AND FROM NEUTRON BETA DECAY. PARAMETERS ARE INVARIANT 

ENERGIES €i, CROSS-SECTION RATIO p, FRACTIONAL ENERGY X, 

AND MULTIPLICITY (. 






XL 



-hi- 



.XL 



-hl- 



1 390 0.65 0.8 1/2 0.05 3/2 0.001 0.5 

2 980 0.44 0.6 1/2 0.05 3 0.001 0.5 

3° 2000 0.23 0.4 1/2 0.05 6 0.001 0.5 

" In statistical multi-pion production at high energies, smaller val- 
ues, e.g., X3 ~ 0.03 can be found via detailed simulations (e.g., 
IMiicke et al.l200Cl : IMurasd200'ff) . with correspondingly higher mul- 
tiplicity, such that the product X3C3 remains at about the same 
value as given here. 



For a continuous jet, the total jet power supplied by the 
black hole is composed of three terms from the magnetic- 
field, particle, and photon power. In the naive standard 
model, we write the total jet power for a two-sided jet as 



2nj/3cR^r^ 
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^EC 



, , , , (37) 

(|Celotti fc FabianllTOQllCelotti fc Ghisellinill2008D . The 
final two terms in this expression represent the photon 
power from synchrotron and Compton processes, respec- 
tively, and are derived in Appendix \^ Here Lsyn and 
Leg are the measured apparent bolometric luminosities 
of the synchrotron and Compton processes, respectively, 
noting that the synchrotron term applies to emission that 
is isotropic in the comoving frame with beaming factor 
(Sq , and so would also apply to synchrotron-self Compton 
emission, while the Compton term applies to the external 
Compton process with beaming facto r 5^ , which applies 
both in the Thomson ([Permed 1 1995(1 and Klein-Nishina 
regimes (jGeorganopoulos et al.ll200ll 120041 ). 

The comoving spherical blob volume is, of course, 
Vl = A-Kr'^ /i, and the particle power, assumed to be 
dominated by the hadrons (in particular, the accelerated 
protons), is just 



£' 



dl'l'N'^ii) 



(38) 



The jet opening solid-angle Q,j = irr'^ / R^ . 

The comoving magnetic field B' required to accelerate 
particles to energy E after escaping the jet is restricted 
by the Hillas (1984) condition, namely that the comoving 
particle Larmor radius r'^ = E'/QB' ^ E/TZeB' < r^. 
The causality condition restricts the blob radius to be 



< 



cSutv 
1 + . 



(39) 



where the measured variability time is denoted by ivar = 
10'^ivar,3 s. Thus the minimum magnetic field B'^-^^ re- 
quired for proton acceleration to energy 10^°i?2o eV is 



given by 

^ {l + z)E ^ E,,{l + z) 

""" recMvar (r/100)(5D/100)ivar,3 ^^ 

The accelerated proton distribution in the fluid frame 
is assumed to take the form 

%^N'^{i,) - K'H{i^ - 7:„iJ7r^ exp(--!^) , (41) 

/max 

so that 7^^^^ ^ 1.1 X 10"^2o/r when B' = B'^^^. The 
normalization for K' in terms of the jet-frame particle 
energy content 



■^par 
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Lpl. 



(42) 

where G = 7,'„ax/7min- The total particle power, as- 
sumed to be dominated by accelerated protons, is from 
Equation ([57]) given by 



/ *,par — „ , pel Cpa^ 



(43) 



The total magnetic-field power V^^^b = ficr'^T'^{B''^ /A). 
Using Equation (gD]) with B' = 5^-^^ gives 



n,B ^^^^8.3 xlO-/3§I ergs- 



(44) 



for relativi stic outflows, which recovers a famil- 
iar result (iWaxmanI 120041: iFarrar fc Gruzinovl 120091: 
iDermer fc Razzaaud I2010D . Note that unlike the 
magnetic-field power, the particle power T'^^par oc r'^^ , 
and becomes larger for smaller <var- 

4.5. Maximum Particle Energy and Bulk Lorentz Factor 
The minimum bulk Lorentz factor from 77 attenua- 
tion is given by Tmin = [crrdKl + zY f^til At^^^m,eC^Y/^ , 
where /e — 10^^^/_i2 erg cm~^ s~^ is the vF^ flux of 
the target photons (assumed to be produced co-spatially 
with the highest energy photon with energy irieC^ti) eval- 
uated at e = 2r2^i„/(l + zfei. For 4C +21.35, 



14.1 



/_i2(^/100 GeV) 

(tvar/600 s) 



1/6 



(45) 



Photons measured with energies of 100 GeV prefer- 
entially interact through 77 processes with photons 
wi th energies E = m,eS? e = 0.5 keV, where /_i2 = 
1 (jTavecchio et al.l 120 111) . This estimate assume co- 
spatiality of the emission, and only takes into account 
77 absorption by photons in the jet. 

The Hillas (1984) condition, with B' given by Equa- 
tion ([3]) and the emission region size scale r[, given by 
Equation (j39p . implies a maximum energy 



E^ 



= Ze^B'r'^. '^ 



4 X 102"Z 



100 



(i + z)6oosy y (Ac/100) 



eV 



(46) 



of particles accelerated by and escaping from the jet. 
Here (5d ~ E is used. Thus, though other losses such as 



photohadronic cooling may be relevant, acceleration of 
UHECR protons is feasible in the inner jet of 4C -1-21.35 
provided F > 100 and B' sa 10 G. There is no conflict 
between this value and Equation (|45|) . 



4.6. Single blob versus continuous jet 

In our formulation of the problem, the secondary pro- 
duction spectrum was derived for a single one-zone blob, 
but the jet power was derived for a continuous jet. The 
energy density in the blob, which is mainly in the form 
of energetic particles, is equated with the product of the 
energy density used to determine the power of the contin- 
uous jet and the blob volume determined through Equa- 
tion ([55)1 . Consequently the single blob represents only a 
single slice with comoving width r[, and stationary frame 
width w ctvar/(l + z) of the continuous jet, and has a co- 
moving particle energy content given by Equation (|43|) . 

This is not, however, the only or the most intuitive 
normalization. As outlined in the estimates. Sections 3.2 
and 3.3, we can alternately constrain the total energy in 
an outburst that lasts for ivar to be Eddington-limited, 
so that the total particle energy £par ^ iEddivar/(l + z), 
or the comoving particle energy 



£par ^ ^Edd^iia 



,/F(l + z). 



(47) 



Comparing with Equation (j43p and requiring 7^*, par < 
Z^Edd shows that the two normalizations differ by a factor 
2(5d/3F, which is of order unity for emission within the 
Doppler beaming cone. 

We first derive the neutron production luminosity for a 
single blob traveling through the external radiation field. 
Because it represents only one zone of the continuous 
jet, it can severely underestimate the possible secondary 
power for a persistent jet. Afterword we consider sec- 
ondary emission from a continuous jet, noting that this 
produces a nonvariable/continuous outfiow with duration 
of ss {l+z)R/c, whereas the single blob gives an outgoing 
pulse of secondaries that preserves the engine variability 
timescale. 

4.6.1. Discrete blob 

The apparent isotropic luminosity of secondary neu- 
tron produced by cosmic-ray protons in a jet traveling 
through a background external isotropic monochromatic 
radiation field characterized by uo/e* is, from Equation 
given by the expression 
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(48) 



restricting to the s — >■ 2 limit in the final expression. 
Here, the term [1 — (yi/e*)^] in Eq. ([SS)) is neglected 
because it makes only a small correction to the integral. 
If a single blob explains the 4C -1-21.35 observations, 
then a minimal requirement is that the apparent sec- 
ondary power in neutrons exceeds the radiant power 



10 



L^ — 10*^^48 erg s^^ measured in the observer direc- 
tion, and which is subsequently converted to an apparent 
isotropic luminosity of ~ 10^^ erg s~^ in 7 rays at the 
pc scale. The condition 47rLjj"<=(f7„) > IQ'^^Lin erg s'^ 
implies, when s ~ 2, 



i^Unc > 



3.6 X IQ^'^L. 
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(49) 



where Ti = 10^^Ti^_i, and A = (2(5De*7max/^i)- From 
Equations (H^ . P5|) and (|i^ . the absolute jet power is 
required to be at least 
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(50) 

The analogous expressions for a thermal graybody ra- 
diation field, from Equation (1311) . is 
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where C = (2e(5D7max/ei)7 implying 
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From Equations p2)) and (l43l) . the absolute jet power 
when s = 2 is at least 



P*. 



par 



> 



3x 10'^'^{l + z)L^ii (r/100)MnG' 

tvar,3(uo/e)4Ti^_l(fc/100)6 1nC 



erg s 



(53) 

Using the alternate normalization for the single blob co- 
moving energy content. Equation (j47p . amounts to mul- 
tiplying these estimates by the 2(5D/3r factor with no im- 
pact on these jet-power estimates for jets observed within 
the Doppler beaming cone. 

Equations ([50)1 and (|53p can be read as saying that to 
produce apparent isotropic neutron luminosity of ~ 10'** 
erg s~^ from an outburst lasting ^ 10'^ s to the observer 
requires a jet with Jd ~ E « 10^ and absolute jet power 
of ^ 10"*^ erg s~* primarily in the form of cosmic-ray 
protons with an s = 2 spectrum, confirming the esti- 
mates of Section 3. The large Doppler factors and asso- 
ciated jet coUimation, with jet opening angle 9 ~ l/E, 
overcome the photohadronic inefficiency and shortens the 
observed timescale of emission to be consistent with sub- 
Eddington luminosities and rapid variability. 

4.6.2. Continuous jet 

The differential secondary neutron production lumi- 
nosity 7„L„(7„,J7„) applies to a single blob traveling 
through a background external radiation field of spa- 
tial extent R. This emission persists for a time At — 
R{1 + z)/(3rS£,c for an observer at angle 6 with re- 
spect to the jet axis, so we can write the time-dependent 
secondary production spectrum for a single zone as 

7nin(7n,^n,0 = lnLn{ln,^n)H{t\tQ,tQ + M) , whcrC 

io is some arbitrarily chosen zero of time. 



For a succession of blobs that comprise the continuous 
jet, the time-averaged spectral luminosity is 

(7n-^n(7n7^n)) = (/0*Ai) 7„L„(7„,fi„) < 

{l + zfR 



civarE^D 



7n-^n(7n,i^n) , 



(54) 



K'TiSl InC , 

(51) 

(52) 



where p* is the rate at which blobs are ejected. The 
largest possible value of p* corresponds to the ejection 
of blobs at the rate /3c/ Ar, where Ar = civar/(l + z) is 
the stationary-frame width of the blob, that is, when the 
time between blob ejection is just Ar//3c, corresponding 
to the continuous jet limit. This explains the coefficient 
in the final term of Equation ([M)) . There is reduction 
by one power in the beaming factor from a discrete blob 
to a continuous jet, as is well-known (|Lind fc BlandfordI 

mil). 

4.6.3. Internal synchrotron photons 

An alternative to assuming large Doppler factors is to 
have small Doppler factors ((5d ^ 100), so that the inter- 
nal sy nchrotron photons present a de nse target radiation 
field (| Atovan fc Dermeii l200ll [20031) : see Appendix p 
But the Doppler amplification of the received luminosity 
would be much reduced by such small Doppler factors, 
so that absolute jet powers larger than the Eddington 
luminosity would be required. Consequently, the large 
Doppler factor solution seems preferable unless an al- 
ternative mechanism to make rapid variability at the pc 
scale is devised. Note furthermore t hat photons from the 
direct accretion-disk radiation field ([Miicke fc Protherod 
120011 ) ■ which provide an additional important target pho- 
ton source that has not been considered here, is more im- 
portant for jetted emitting regions far from the black hole 
when the Doppler factor is larger , because aberration of 
accre tion-disk photons is greater ([Dermer fc Schlickeiseii 
[2OOI . 

4.7. Calculations 

We use the formalism developed here to make calcu- 
lations of secondary neutron and neutrino production 
spectra from the interactions between a power-law dis- 
tribution of UHECR protons, isotropic in the jet rest 
frame, that interact with BLR radiation in the inner 
jet to make secondary neutrons, pions, and neutrinos. 
The results are expressed in the form 4:TrEL{E, ft) — 
4:'KE(PL{E,n.)/dEdn. (units of erg s~*), so as an ap- 
parent isotropic luminosity for the production of secon- 
daries with that energy. Given that we are scaling to 
4C -h21.35 with 47rd| = 6.52 x 10^^ cm^ and bolometric 
energy fiuxes during fiares exceeding 10~^ erg cm~^ s~*, 
this corresponds to apparent bolometric 7-ray luminosi- 
ties > 10'*® erg s^*. 

Flaring results to model the SED of 4C -^21. 35 from 
photohadronic neutron and neutrino emissions made by 
UHECR particle acceleration in the inner jet are shown 
in Figures [3] and [H respectively. In the picture studied 
here, cosmic-ray protons in a single-zone blob with F — 
100 and (5d — 200 — so looking directly down the jet — 
interact with photons of the BLR and the dust torus. 
Here we normalize to a total jet power T'* = 2 x 10*^ erg 
s~*, which is the Eddington luminosity of a 1.5 x IO^Mq 
black hole. We use ivar = 600 s for the variability time 
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E (eV) 

n 

Fig. 3. — Black solid curve shows the secondary production 
spectrum of neutrons made from cosmic-ray protons accelerated in 
the inner jet of a blazar. The jet is assumed to have bulk Lorentz 
factor r = 100, be aligned along the jet axis {Oohs = 0°, 5d = 
200), and interact with BLR and torus photons. The separate 
components from interactions with Lyo, scattered accretion disk, 
and warm and cool dust radiations are labeled. The contributions 
from the different resonance and multipion interactions are shown 
separately for the scattered accretion-disk radiation component. 



scale, and Equation (H71) to normalize the energy content 
of jet protons. This jet energy is fed into magnetic field, 
photons, and particles, limited by V^, = L^dd- To ac- 
celerate protons to 10^" eV, B' > 1 G, from Equation 
(|40|). Most of the power goes into particles; the mag- 
netic field required to accelerate 10^° eV protons leads 
to a total jet magnetic-field power of « 10^^i<^|o ^rg s~^, 
from Equation (|34|) . which is « 5% of iEdd- The photon 
power is smaller still, due to the small beaming factor 
w 1/r^. We let 7'niin = r, corresponding to the mini- 
mum Lorentz factor of protons swept into the jet. The 
effective photon densities for the Ly a and scattered ac- 
cretion disk radiation field are given by Equations ([5]) and 
([6l) , respectively. The total power that is not in magnetic 
field or photons is assumed to be mostly in the form of a 
nonthermal particle distribution of the form of Equation 
dU), with s = 2.1. 

Figure [3] shows the neutron production spectrum for a 
head-on F = 100 jet, 6'obs = 0°. The large apparent neu- 



tron luminosities, exceeding 10 erg s~ from an AGN 
with absolute jet power of just 2 x 10**^ ergs, shows that 
photohadronic pr ocesses in the inner jet can be very effi- 
cient (contrary to iNalewaiko et al.ll2012D , with the warm 
dust more effective than accretion-disk radiation for hard 
cosmic-ray proton spectrum that extends to 10^° eV. The 
two peaks found in the production/injection of neutrons 
into the pc-scale region where only IR dust photons are 
present means that only the highest energy neutrons with 
En > 10^® eV, are effective at making photopion secon- 
daries beyond the inner jet within the BLR. 

Figure U] is a calculation of the neutrino production 
spectra for an inner-jet system with the same parame- 
ters as in Fig. [3] No mixing corrections are made. The 
apparent bolometric luminosity is ^ 10^^ erg s~^, which 
as noted before, requires large Lorentz factors and pref- 
erential Doppler factors, and production in the inner jet 
for short variability. The components labeled by tt — > z/ 
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Fig. 4. — Black solid curve shows the secondary production 
spectrum of neutrinos made from cosmic-ray protons accelerated 
in the inner jet of a blazar, using the same parameters as in Fig. 
[3] The separate components for neutrino production from inter- 
actions with Lyo, scattered accretion disk, and warm and cool 
dust radiations are labeled for the channel (p -I- 7 — >)7r — >■ u. The 
corresponding line styles are shown for the production of /3-decay 
neutrinos denoted hy n —^ u, which is formed when the neutrons, 
assumed to escape the production site and travel rectilinearly until 
decay, make neutrinos through neutron /3 decay n — >p-|-e~ +Ve. 



refer to neutrino production in photopion interactions, 
p -|- 7 — ^ TT — >■ J/. The low-energy feature in the neutrino 
production spectra labeled n — >■ 1/ are /3-decay electron 
neutrinos, when the bulk of the neutrons escape and de- 
cay in intergalactic space. Both the neutron and neutrino 
spectra become very steep, harder than +2 in EL{E), 
below PeV energies, as a result of thre shold effects. De- 
tailed calculations (e.g., iMurasd '20071 find softer spec- 
tra below PeV energies than calculated here due to the 
spread in the energy distribution of low-energy pions, es- 
pecially from 7r//i decay kinematics and partially from 
multi-pion production. 

Fig. [5] compares the apparent isotropic luminosities 
(or production spectra) of neutrons (blue curves) with 
neutrinos (red curves). Using the same parameters as 
Fig. [21 this figure shows production spectra calculated 
for F = 100, 50 and 20 at observing angles 9 = and 
9 = 1/F. The comparison reveals the strong, S^ beam- 
ing factor of this process (Eq. I35p when viewing off-axis, 
leading to a factor 2^ = 32 reduction of the production 
luminosity at the Doppler beaming angle compared to 
the luminosity along the jet axis. The relative fluxes as 
a function of F at the same observing angle (which would 
be peculiar except for = 0) is closer to ex F^, due to 
kinematic effects and normalization. Neutrons are made 
as secondaries of protons with En ~ O.SEp and carry the 
bulk of the energy of the incident protons. Neutrinos, 
by comparison, carry only « 5% and /3-neutrinos only 
^ 0.1% of the energy of the original neutron. 

Figs. [3] ~ m apply to single one-zone blobs traveling 
through a beam production region of extent i?BLR = 
2 X 10^ cm. The duration of the neutron production 
for on-axis observers is ~ (1-1- z)i?BLR/2F^c ~ 500 s, 
and shorter if the accelerated proton distribution decays 
before it traverses the length of the jet. 

Calculations of secondary neutron and neutrino pro- 
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E = E (eV) 

Fig. 5. — Comparison of production spectra of neutrons (blue 
curves) and neutrinos (red curves) for jets with different Lorentz 
factors and observing angles, as labeled. For the standard inner 
jet and accelerated cosmic-ray spectrum described in the text, the 
production spectra are calculated for F = 100, 50 and 20 at obser- 
vation angles 6 = and d = 1/F. 




s =ee or v 

1 



-1® 



Fig. 6. — The function (p(s)/s^ , calculated from Equation I ISSH . is 
given by the solid red curve, along with asymptotes in the low- and 
high-energy limits from E qua tion l|59| l. The dark blue curve is the 
func tion T(u), Equation II63I I. with asymptotes given by Equation 
ll64t . 



duction spectra from a continuous jet allow for larger 
secondary luminosities at the expense of variability. On 
the other hand, smaller values of F might be more typical 
during these extended periods, and it is less likely that 
we would be looking directly along the jet axis. Even for 
smaller jet Lorentz factors F sa 20, detection of neutron- 
induced or neutrino emission might be compensated by 
its longer duration, though compromised by the larger 
backgrounds for neutrino detection for the longer time 
windows. 

5. PAIR-PRODUCTION OPACITY AND PHOTOPION 
EFFICIENCY FOR ISOTROPIC RADIATION FIELDS 

Our results to now are fairly intuitive. Cosmic-ray pro- 
tons with energies 7p > ethr/e* interacting with an ex- 
ternal radiation field with mean photon energy nieC^c* 
have a photohadronic neutron-production efficiency 



1 - 10% in the inner jet consist- 



Kp^apjRuo/meC ei, 
ing of the BLR and IR dust torus radiation field local to 
the supermassive black hole in 4C -1-21.35 (equation (fTUI) 
and Sections 3.2 and 3.3). With the stipulation that the 
absolute jet luminosity is Eddington-limitcd and that the 
bulk of this power is transformed to an s = 2 spectrum 
of protons with maximum escaping energies ~ 10^^ eV, 
then an impulsive jet has to have (5d ^ 100 to produce 
a highly collimated beam of outflowing UHE neutrons 
with apparent luminosity sa 10^^ erg s~^ in a pulse of 
particles < lO^'^ cm in width measured in the black-hole 
frame. The large Doppler factor coUimates the beam and 
preserves the rapid variability. For a continuous jet, ap- 
parent isotropic neutron luminosities of sa 10"*^ erg s~^ 
are formed by a jet with more modest Doppler factors. 
The outflowing neutron or neutrino fluxes will persist for 
longer times, though without short timescale variability. 
We now turn our attention to the 77 opacity and pho- 
topion cfhciency of 7 rays and neutrons in the radiation 
field formed by the BLR and dust torus, for application 
to 4C +21.35 or other blazars. We work under t he as - 
sumption of local isotropy; the paper by iGouldl (|19790 
calculates energy density with this assumption relaxed. 



5.1. 77 Opacity 

The opacity T^j{ei) for a 7 ray with energy meC^ei to 
travel a distance R through an isotropic radiation field 
with energy density uq = nieC^ L deenph(e) consisting 
of photons with energy meC^e described by the number 
density distribution riph(e) is given through the relation 



dr. 



77 



dx 



^2 



de 



^1 Jlje 



'^ph(e) V'(so) , (55) 



where sq = eei, 



/So 
ds 



sayj{s) 



„2 ' 



(56) 



and aj~f{s) is the 77 pair-production cross section. Thus 
Tir^R 



T^^Ul^ 



Jl/ei 



dee '^nph{e)ip{so) 



(57) 



The function 

(p(so) = (2so -2 + —) In Wo + 2(1 - 2so)\/l-So^ + 



-. 00 

\nwo[4:\n{l+wa)-3\nwo]- -TT^ +4:Y^{-y'-^n~^w^ 



n=l 



'0 ' 



(58) 



where wq = 2sq\Ji — Sg ^ — 1. The asymptotes of '■p{sq) 
are 



(2so-Mn4so)(ln4so-2)- 



7r^-9 



(p{so) -^ 



-hSo^(ln4so + 9/8) + ... so > 1 



|(S0- 1)3/2 + |(.„-1)V2_ 

(253/70)(5o- 1)^/2+ _, 



So - K 1 

, , , (59) 

(iGould fc Schreder I 119671: iBrown. Mikaelian. fc Gouldl 



13 



The pair-production opacity for a photon with energy 
TOgC^ei to pass through an isotropic monochromatic radi- 
ation field with target photon energy TTigC^e* is therefore 



r(ei) 



v{so) 



where r^ 



tttIRuq 3 gtRuq 



ruec^ei, 



e^nieC^ 



(60) 
A graph of the function ip{sq)/sq used to calculate r^^i^ci) 
in an isotropic photon field, including asymptotes from 
Equation (|59)) at small and large values of sq , is shown in 
Fig.[6l Note the ln(ei)/ei dependence at so ^ 1 because 
'^(■5o)/so ~^ 21n(0.54so)/so in the limit sq > 1. The 
function (p(sq)/sq reaches a maxiumum value of « 0.56 
at So — 3.54, and falls by a factor of « 20 when sq 
goes from 3.5 to 300. For Ly a line photons, the peak 
opacity corresponds to about 90 GeV in the stationary 
frame. Neglecting photon broadening processes, there is 
a sharp lower limit to the opacity at ei = 1/e*, that is, 
at 25.6 GeV in the stationary frame for opacity from Ly 
a photons. 

Substituting the graybody spectrum, Equation (|29p. 
into Equation (1551) gives 



r|!;(ei) 



2a)gR 
Ac 



2algR „ 



de 



Viso) 



Q^F{v) 



i/ei exp(e/e) - 1 
15r^R uq 



TT"' JTleC 



•e 



F{v) (61) 



for the energy-dependent 77 opacity through a graybody 
radiation field, where the fine-structure constant a/ = 
e'^/hc ^ 1/137, and 



1 
7^ 



(62) 



Here g is the graybody factor representing the ratio 
uo/uhh of energy densities of the graybody and black- 
body radiation fields with temperature T (see Eq. 
and the blackbody 77 opacity function 



Tiv) = v' 



dw 



exp(w) — 1 



(63) 



This function reaches its maximum value at J-pk = 
F{vpk) = 1.076 at Vpk = 0.51. The asymptotes of J-ii^) 
are 

r V^ exp(-j/)(l -K I;) , i/>lorei<l/e 

[ ^ ln(0.47/i/), i^ < 1 or £1 > 1/9 

(64) 
Fig. [6] is a graph of the function J^{v) along with its 
asymptotes. 

5.2. Photopion Efficiency 

The derivation of the photopion energy-loss rate of an 
ultrarelativistic neutron with energy En = 7nm„c^ = 
10^°i?2o eV in an external isotropic radiation field de- 
scribed by the function npii(e) follows closely the deriva- 
tion sketched in Section 3.3 for cosmic-ray protons in a 
jet, and is given by 



*07(7«) = ^ 



de 



^ph(e) 



27„e 



d€r era^j{er)K^j{er) 




CO 

Fig. 7. — The function I (to) and asymptotes, from Equation I I36II . 
that describe the photopion energy loss rate of protons or neutrons 
with photons of a blackbody or graybody radiation field. 
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de nph(e) 
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27«e 



(65) 



using the approximation cr0-y(e,.)i4r^-y(e,.) = aH{er — lt\i-c) 
where, as before, ethr — 400 is the photopion threshold 
photon energy and ct = 70 ^b is the product of the pho- 
topion energy-loss proto n cross section and inelasticity 
(| Atovan &: Dermerl I2003D . In a graybody photon field 
given by Equation ([M]) with graybody factor g, 



tA^iln) 



SiTcgcrQ^ 



^H, 



(66) 



where /(w) 
parameter 



12(0;) -uj^Io{uj), from Equation ^^. The 



ethr 

2^ 



(67) 



characterizes the different regimes of photopion interac- 
tion of an ultrarelativistic cosmic-ray proton in a black- 
body photon distribution, and whether the proton inter- 
acts with the exponentially declining number of photons 
in the Wien regime (w > 1;7„ < 200/8 or En < 10^^ 
eV/(rdust/1200 K)), or with the bulk of the distribution 
(uj « 1;7„ » 200/8; ^„» lO^^ eV/(Tdust/1200 K)). 
Expressing the photopion efficiency 



V<pj — ^Ic/i^' 



■7 5 



where the light-crossing time scale tic 
region of size R, then 



(68) 
R/c across a 



.gb 



vi;iEn) 






lico) 



15 aRuo ^, , 



TT^ nieC 



29 



(69) 



The photopion efficiency for monochromatic line radia- 
tion is simply 



v't^'iEn) 



Uq&R 



{1-X-^)H{X-1) , X: 



27Tie* 

ethr 



(70) 

and approaches a constant value at 7„ 3> ethr/'2e* (com- 
pare Equation (fTH)) ). We show jjpj in Figure 9 using 
^BLR = 2 X 10^^ cm. One should keep in mind that 
i? = 1 pc should be used for the photomeson production 



14 



between the neutron beam escaping from the BLR region 
and target photon fields from the dust torus, and beamed 
neutrons are essentially depleted in the dust torus. 

5.3. Ratios 

For isotropic monochromatic radiation, the ratio of the 
peak 77 opacity to the photopion efficiency in the limit 
of high energy neutrons is 



n 



line 



pk^ 



A^T) 



V<t>-^{x > 1) 



0.56 



7rr| 
a 

^2000 



v{sq) 



pk 



(71) 



The ratio of neutron energy at a; = 1 and photon energy 
at the opacity peak is 



e 



line 
P7 



Enjx = 1) 
ppk 



Cthr 



10^ 



£;f"'' me (2 X 3.54) 

For thermal blackbody or graybody radiation, 

7rr|^.,„„„„ 1.076 
CT Ipk 



n 



thermal 



3600 



r(3)C(3) 



1600 



(72) 



(73) 



and the ratio of neutron energy at 2; = 1 and photon 
energy at the opacity peak for the thermal photon dis- 
tribution is 



c, 



thermal 
P7 



4771, 



■ ethr ~ 2 X 10^ 



(74) 



Corresp onding relations in relativistic jet sources are 
given bv lDermer et al.l (120071 ). 




E^ = (1+z) E(eV), E^(eV) 

Fig. 8. — Pair production opacities of 7 rays and photopion ef- 
ficiencies of UHECR protons or neutrons in a monochromatic Ly 
a radiation field in the BLR, with the radiation field scattered by 
thermal electrons in the BLR, and with the dust radiation fields of 
the IR radiation field found both in the BLR and the dusty torus. 
Photon energies E{cV) and neutron energies _E„(eV) are measured 
in the local source frame. For 7 rays, R = Rblr = 2 x lO^'^ cm is 
used for the emission radius of the BLR region and /? = 1 pc is used 
for the emission radius of the dust torus. tjHE neutron production 
is assumed to occur in the BLR region, where both BLR and IR ra- 
diation fields are important, and escaping neutrons are assumed to 
interact both with Ly a and scattered radiation within R = -RblRi 
but only with the IR radiation of the torus within R = 1 pc. Light 
dotted curves show opacities and photopion efficiencies using the 
graybody approximation for the scattered accretion-disk radiation 
field described by Equation l(6)l . 



5.4. Application to 4C +21.35 
From eqs. ([5]) and (|60|. we find that the quasar BLR 



opacity to Ly a photons is 



T^T,(ei) ^ 4Oi?i70_i 



v{so) 



(75) 



The Ly a opacity with i?i70_i = 2 is shown in Figure. [51 
The large opacity in the VHE range precludes the pho- 
tons detected with MAGIC from 4C +21.35 from being 

made in the BLR. 

Using the iMalmrose et al.l (|2011l ) parameters for the 
quasi-thermal IR emission from the hot dust, the gray- 
body factor 






0.18 



Lie I 



i?p,2 VI2OOK 



(76) 



taking uyr ^ L/AtiB?c, and noting that the energy den- 
sity Ubb(1200 K) = 0.0156 erg cm^^ of a 1200 K black- 
body radiation field. From Equation ((6T|) . 



bb 

77 



(ei) 



2a)gR 
Ac 



Q^F{u) 



2000 ^F{v) , 

^pe 

pk 



(77) 



with a peak opacity at v ^ 0.5, implying E\ = 5 TeV 
For the warm dust component at T = 660 K, 



^bb 



(ei) 



350 ^J-(:.), 

^De 



(78) 



with a peak opacity at Ef~ = 10 TeV. 

In addition, there is the scattered radiation field from 
free electrons in the BLR that scatter radiation emitted 
by the optically thick Shakura-Sunyaev accretion disk 
and a hypothetical hot X-ray pair plasma found close to 
the black hole. The optical depths to these photon fields 
are calculated using Equation ((57)) and the relation 



„c^e^nph(e) = eusc(e) 



eLie) 



(79) 



AttR^c 

for the scattered radiation field. The Shakura-Sunyaev 
disk spectrum in the form 

-^disk 



e^disk(e) 



r(4/3) 



4/3 



exp(-e/e„ 



(80) 



is used to estimate T^-y(ei) from its scattered radiation 
field. For a hot plasma making an X-ray spectrum with 
index (a = 1/2) and luminosity Lx, 



eL-(e) 



Lx 

r(i/2) 



1/2 



exp(-e/emax) i?(e- emin) 



The normalization becomes increasingly accurate when 
Fig. [S] shows calculations for a disk spectrum with lu- 



minosity Ldisk = 2 X 10"*^ erg s ^ and i?n.„.., 

eV, and an X-ray corona with Lx — 4.5 x lO**^ 



50 
erg 



15 



s~^ and iJ^max.JC — 5 keV, us ing a dust covering factor 
of 22% (jMalmrose et al.ll2011[ ) and taking the scattering 
optical depth Tsc = 0.01. The photopion efficiency was 
determined for the scattered radiation fields by substi- 
tuting Equations (|80l) and (ISTI) into Equation (|55|) . using 
Equation ([7^ and Equation ([BS]). For comparison, the 
dotted curves show the opacity and photopion efficiency 
when the scattered accretion-disk radiation field is ap- 
proximated by a thermal graybody spectrum using the 
parameters from Equation ([6]). Note that radiation fields 
with higher effective temperatures but comparable lumi- 
nosities provide negligible additional opacity because of 
the smaller density of target photons. 



6. 



SYNCHROTRON RADIATION FROM UHECR 
SECONDARIES 



6.1. Analytical considerations 

The results in Figure [8] show that the multi-GeV 7 
rays detected from 4C -1-21.35 with MAGIC cannot be 
made in the BLR, and that inner jet radiation should be 
strongly attenuated above several GeV from scattered 
accretion-disk radiation. For an accretion-disk power of 
5 X 10^^ erg s~^, Rn sa 2 from Equation ([1]), and if 
the accretion-dis k power is even larger, as indicated by 
the Swift data ([Tavecchio et al.l 1201 ID . the 77 attenu- 
ation is likely to be even larger than shown in Fig. [8l 
Thus the emission region for VHE photons must at least 
be situated beyond the edge of the BLR, but even then 
photons with energies > 200 GeV are strongly absorbed 
by the IR radiation from the torus. The actual attenua- 
tion depends precisely on deviations of the dust spectrum 
from a graybody in the Wien regime, but even with the 
exponentially declining number of photons in the near 
infrared spectrum, there is still large opacity for VHE 
emission above a few hundred GeV. The lack of attenu- 
ation in the MAGIC data to - 400 GeV (~ 600 GeV in 
the source frame) implies that the VHE 7 rays have to 
be made at and beyond the pc scale. 

Escaping neutrons easily leave the BLR but suffer 
strong photopion losses when passing through the in- 
frared field of the dust torus, where neutrons would 
deposit a large fraction of their energy during escape. 
Hence, we should expect formation of beams consisting of 
UHE 7 rays, pairs and neutrinos. The photopion produc- 
tion efficiency riff,^(jn) of an escaping neutron or proton 
with Lorentz factor 7„ passing through, in particular, the 
hot dust radiation field of 4C -1-21.35, is estimated to be 
V'^-riln) ^ 2.5(uo/e*)4-Rpc (see Equation ^0^ and Fig.[8l 
noting (e*) = 2.70), and this energy is transformed into 
ultra-relativistic e^-e^ pairs, photons and neutrinos. 

In the case of 4C +21.35, efficient photopion losses 
apply to neutrons or protons with w ^ 1, which for hot 
dust with Q = 2 X 10"'' implies a range of energies 
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= 6 X 10^^ eV < E„ < Em.^ = 



IQ^^E2Q eV 



(82) 

for which photopion losses are large, as can be seen from 
Figure [51 

In each n-l-7 — )■ p-|-7r~ interaction, roughly 20% of the 
neutron's original energy is given to two pionic photons, 
« 5% is given to an electron, and ~ 15% is given to three 
neutrinos. Secondary electron energies therefore range in 
value from sa 3 x 10^^ eV to 2 x 10^^i?2o eV, implying 




Fig. 9. — Energy-loss mean-free-path (MFP) for synchrotron and 
Compton losses of an electron with Lorentz factor 7e. Synchrotron 
loss MFPs are calculated in a tangled magnetic field ranging from 
1 /iC to 10 mC, as labeled. Compton energy-loss MFPs are plotted 
for the CMB radiation, and for graybody IR radiation fields from 
hot (T = 1200 K) and warm (T = 660 K) dust, with a graybody 
fac tor correspondin g to a distance il = 1 pc and properties given 
by IMalmrose et al] l|20111) . Vertical lines represent the range of 
injection electron Lorentz factors for electrons formed by the decay 
of 7r+ from photopion interactions of UHECR neutrons. The cross- 
hatched region represents the range of 7e and magnetic fields where 
10 minute variability can be preserved (see text). 



Lorentz factors in the range 



6 X 10^° < 7e 



(83) 



so 0.6 < 711 < 100£'2Q. The synchrotron photons from 
these hyper-relativistiqj be produced through secondary 
e+ and e~ have observed energies 



-^syn — TTLqC 6g 



Smec'^Bj'^ 

2Bcr(l + Z 



120S^G,pc7ii MeV 



(84) 

where B^r = 4.414 x lO^'^ G is the critical magnetic field 
and Bpc is the magnetic field at the pc scale. Thus the 
photon energies of the first generation synchrotron spec- 
trum are in the range 



43 MeV < 



E. 



'syn 



B 



7j.G,pc 



< 1.2E^o TcV 



(85) 



The second-generation synchrotron spectrum from e+-e~ 
pairs made when synchrotron 7 rays form pairs through 
77 attenuation is important for larger, ~ niG fields. If 
a synchrotron 7 ray forms a pair, with each electron 
and positron receiving half the energyQ then the second- 

^ The se electrons are "hyper-relativistic" in the sense 
(Dermer .STAtoyan 2004) that such high-energy electrons cannot 
be directly accelerated by traditional Fermi acceleration mecha- 
nisms that compete directly with synchrotron losses. 

^ This is not a good approximation when the photons pair- 
produce in the Klein-Nishina (KN) po rtion of the pair-p roduction 
cross section (see, e.g., Murasc 2009; Kotera et aI.|[20Tll 'l. Escape 
of UHE 7 rays depends on the internal target photon spectrum, 
which may be suppressed by self-absorption of the synchrotron 
spectrum or the Rayl eigh- Jeans portion of the black-lsody spec- 
trum (Razzaque et al.ll2004 : Murase 2009) . 
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generation synchrotron spectrum ranges between 



0.052 eV < 



-^syn 
-°mG,pc 



< 35EL MeV 



(86) 



However, those second generation synchrotron photons 
are typicahy irrelevant for the highly variable emission 
(see below). 

Synchrotron losses will dominate when the synchrotron 
energy-loss MFP is shorter than other energy-loss MFPs. 
Fig. |9] shows the synchrotron energy-loss MFP in com- 
parison with the Compton energy-loss MFP on the dust 
radiation fields and, for comparison, the CMBR field at 
z = 0.432. The electron synchrotron energy-loss rate as- 
sumes a randomly oriented magnetic field on the pc scale 
of strength i?pG,pc- The Compton energy-loss MFP of 
an electron passing through a graybody radiation field is 
approximated in Fig. [9] by the expression 



Ac (7) = c 



7c 
7 



45A? 



106fl„ 



327rSg(TTe*7 ~ L 
1. 52x10-=' fl^^Ti2oo 



2Aa7 



■pc, 



pc. 



7<7 



7 < 7 < 7* 



7.4x10'' 



_ 4.92 



7r3gcrTe2 1n(0.557e) ~ 
6.27iiTf,ooflg, 
>^ L46ln(l.lxlO''7iiTi2oo) ^^' 



7 >7* = 



2.4x10' 
^"1200 



where the energy density of the blackbody radiation 
field with temperature Q is given by Equation ([50]) . and 
T1200 = T/1200 K. This simple expression connects the 
Thomson and extreme KN asymptotes with a constant 
MFP in the intermediate regime determined by equat- 
ing the Thomson MFP with the value of the KN MFP 
evaluated at the electron Lorentz factor 7* given by 
ln(0. 557*6) = 1. 

From Figure [HI one sees that synchrotron losses of 
electrons formed by neutron photopion production dom- 
inate Compton losses when Bpc > 3/iG, and Equa- 
tion (1851) shows that first-generation synchrotron pho- 
tons with energies exceeding a few GeV are formed when 
^pc ^ lO^G. We now show that such a process can 
accommodate the short variability timescale observed 
with MAGIC if the inner engine generates pulses of 
UHE neutrons that are modulated by activity of the in- 
ner jet on short timescales. In this case, synchrotron 
production at distance R from the central engine can 
vary on timescale At if the emission is made within the 
angle 9t of the line of sight according to the relation 
(1 -|- z)R{l — cos9t)/c — At. This is satisfied for elec- 
trons produced within the angle 



7dfl 



<f).= 



I 2cAt 
R{l + z) 



2.9 X 10" 



' (At/600 s) 



il + z)R 



pc 



The electrons making the GeV/TeV synchrotron radia- 
tion are deflected by the angle ^dfl — \/2/3Asyn/rL = 

6-\/2/3TTe/(JTBj^, where Agyn = QTrmeC^/cTTB^J is the 
energy-loss length for synchrotron losses and tl = 
TrieC^^/eB is the electron Larmor radius. The criterion 



that e'dfl < 9t (|Murasdl20il ) impfies 



B^.G..pcln > 330. 



' (At/600 s) 



Rr, 



(89) 



to preserve observed variability on a timescale shorter 
than w 600 s for 4C -)-21.35. Together with Equations 
(|55)) and (|55)) . this variable synchrotron radiation will be 
detected at photon energies 



4Q /^At^OM < ^^y^^(GcV) < 1.2 X lO^S^o^MG^pc • 

(90) 
The constraints will be relaxed if the coherence length 
of the magnetic field is small in comparison with the 
synchrotron cooling length. 

A number of important points can be made regard- 
ing this result. The synchrotron radiation frequency de- 
pends only the product Bj^ and not on i? or 7 sepa- 
rately. Hence, At essentially depe nds on ii/svn, w ithout 
the explicit dependence on B and 7 (|Murasell2012[ ) . Thus 
highly variable synchrotron emission is expected from 
VHE synchrotron 7 rays. Because the photohadronic 
production kinematics limits the lower bound of the sec- 
ondary electron distribution to be at 711 ~ 0.6, from 
Equation (j83| and Figure [U the characteristic magnetic 
field at the pc scale is implied to be < 3 mC, other- 
wise the energy range of synchrotron photons is out of 
the MAGIC band. Lower magnetic fields are allowed as 
long as the synchrotron cooling length is short enough, 
where lower B corresponds to higher 7 at given -Esyn- A 
definite prediction of this model is that the variability 
times of the 7-ray synchrotron radiation become longer 
at lower energies, because a larger angular range of cool- 
ing and defiecting electrons can then contribute to the 
observed emission. Emission from the inner jet can, how- 
ever, dominate at GcV energies, so this prediction applies 
only to the VHE synchrotron 7 rays. Inverse-Compton 
cascade radiation should be more slowly variable because 
Air > Asyn, and second-generation synchrotron emission 
(Equation ((86|) ') does not contribute to the highly vari- 
able emission since it is typically expected at lower ener- 
gies. 

6.2. Numerical Results 

To demonstrate the resulting secondary spectra pro- 
duced by the neutron beam launched from the BLR re- 
gion, we also perform numerical calculations. In order 
to evaluate photon and pair yields from the photome- 
son production, we use SOPHIA ()Mucke et all I2000D 
and solve t he kinetic equations for injected photons and 
pairs, as in IMurase et al.l (j2012[) . We calculate cascades 
taking into account synchrotron and inverse-Compton 
emissions, and we focus on beamed emissions such that 
^dfl < Ot to evaluate variable emission components pro- 
duced in the dust torus. In this work, we only show the 
highly variable emission component with At < 600 s, 
and we do not include emissions with longer timescales, 
though they, as well as radiation from the in ner jet, 
parti ally contribute to the overall received flux (jMurasd 
[2011) . 

Figure [TOj show results for the case where the neutron 
beam is produced by photomeson interactions between 



17 



erg s 



UHE protons and external photons in the BLR. The 
spectrum used here is the same as that shown in the 
inner jet calculation of Figure [3l which includes target 
Ly a, scattered accretion-disk, and cool and warm torus 
radiation as target radiation fields. Then, we calculate 
the photomeson production by UHECR neutrons in the 
external photon fields of the dust torus, whose radius is 
set to i?dust = 1 pc. The magnetic field is set to 10 ^G. In 
the calculations, the dust temperatures are set, as before, 
to Tdusti = 1200 K and Tdust2 = 660 K, with luminosi- 
ties Ldusti = 8 X 10''^ erg 8"^ and Ldust2 = 10"^^ 
respectively!! 

We found that the required absolute jet power is 
~ 10^® erg s~^, which is consistent with analytical ex- 
pectations. From Figure m one expects that UHE pho- 
tons might escape from their emission regions, where 
they cascade in larger-scale magnetized regions, includ- 
ing elliptic al galaxies, ga lactic winds, galaxy clusters and 
filaments (jMurasd 120121 ) . Since strong radio emission 
is presumably produced by kpc jets in FSRQs, UHE 
7 rays can also be quickly depleted by such an addi- 
tional target radio field, whose scale and intensity is 
uncertain. For demonstration, we therefore consider 
an additional radio field which wou l d be provided by 
the kpc jet. From iTavecchio et al.l (|2010f) . the radio 
field is optimistically assumed to have a typical size of 
^radio = 10^^ Cm and a characteristic vL^ radio lumi- 



nosity = lO'^'' erg s"^ {i^/W Hz) at i^ < W^ Hz, with a 
naive extrapolation to lower frequencies. 

To calculate this contribution, we also evaluate cas- 
cades initiated by UHE photons leaving the dust torus, 
using cAi « (1 -f z)(l — cos6'dfl)min[A-y^,-Riadio]- Rapidly 
variable emission is still possible for the UHE photon- 
induced emission as long as A^-y is short enough ( Murasa 
[2012). The steady VHE cascade emission induced by the 
kpc scale jet is predicted, from Fig. [101 to be at a lower 
flux level than the rapidly variable VHE emission formed 
at the pc scale. The search for the low-level, high-energy 
plateau emissi on is feasible with the Cherenkov Telescope 
Array (CTA) (jActis et al.ll2011[ ) in development. Besides 
the UHE photons and leptons generated by neutron pho- 
tomeson processes at the pc scale (thin solid curves) , the 
subsequent cascade 7-ray spectrum formed in the dust 
torus is shown by the dashed curve. The UHE photons 
escaping into the kpc jet can experience further cascades, 
which makes an additional high-energy component with 
rapid variability. 

The results are not very sensitive to B as long as the 
magnetic field is in the range 10 /iG ^ B < mG, as 
shown in Figure [TT] Just for comparison, in Figure [T^ we 
also show the case where the neutron beam is produced 
by photomeson production between UHE protons and 
nonthermal photons that can be generated by electrons 
accelerated at internal shocks. Here, we inject beamed 
neutrons with a power-law spectrum, Nn{En) ex E~^-^. 
Such a hard spectrum can be realized if the target photon 
spectrum in inner jets has /3 = 2.2 and the proton spec- 
trum has s — 2.0 (App. B). The target photon spectrum 
is taken not to contradict the observations. As in Fig. 
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Fig. 10. — 7-ray spectra formed by a neutron beam generated in 
the inner jet, based on the calculations shown in Figure[3] MAGIC 
data are also overlaid, which have been dcabsorbed using a low- 
intensity model of the EBL (.Dominguez et al..,2011.) . 
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Fig. 11. — 7-ray spectra formed by a neutron beam generated 
in the inner jet, based on the calculations shown in Figure |3] but 
for different values of B. The required absolute jet powers are 
9.5 X 10-*5 erg s"!, 6.8 x lO-^^ erg s"! and 4.1 X 10* erg s"! (top 
to bottom). 
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Fig. 12. — 7-ray spectra formed by a neutron beam generated 
in the inner jet, based on a power-law neutron spectrum with 
s = 0.8 at IQls eV < En < lO^^^'^ eV. We use J dEn Ln{E„) = 
10*'' erg s~^. 



^ For numerical calculations, we use TsciJdisk = 5 X 10** erg s ^ 
and TscLx = 4.5 X 10** erg s'^, respectively. The EBL, CMB 
and cosmic radio background are also included, though they are 
irrelevant for source cascades. 



[TTl the dashed curve represents the contribution of cas- 
cade gamma rays developed in the dust torus, while the 
dotted curve is for the additional kpc jet component. In 
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Fig. 13. — Neutrino and cascade 7-ray spectra formed by a neu- 
tron beam generated in the inner jet, based on the calculations 
shown in Figure |3] Thick lines are for 6(jfl < St, while thin lines 
are for d^fi < 9j = 0.01. Solid curves show the cascade 7-ray 
spectrum formed in the dust torus, where the kpc jet component 
is not included. Note that 7 rays with 9^fi > dt have longer time 
scales, and the flux is suppressed at lower energies. 



this case, however, considerably larger absolute powers 
are required in order to explain the observations. 

In Figure 1131 neutrino spectra are shown. Without 
performing a detailed calculation of neutrino detectabil- 
ity using the IceCube detector response, estimates show 
that a fluence of « 10~^ - 10~^ erg cm~^ in 7 rays as- 
sociated with photopion neutrino production at PeV en- 
ergies is requi red for detection of a fe w PeV neutrinos 
with IceCube (jAtovan fc Dermer|[2001[) . Since the differ- 
ential high-energy neutrino flux of ~ 10~^ erg cm~^ s~^ 
is associated with photopion losses, then flaring episodes 
of days to weeks would be required to detect neutrinos 
with a km-scale neutrino telescope. Although the ex- 
pected fluence is not enough to expect neutrino detec- 
tion by IceCube during a single flaring event from 4C 
-1-21.35, stacked flaring episodes of many such events, or 
searches for PeV neutrinos from 7-ray bright FSRQs like 
3C 454.3, 3C 279, or PKS 1510-089 with high-energy 
neutrino telescopes are crucial for testing the model. 



7. DISCUSSION AND SUMMARY 

The sources of the UHECRs are unknown, but ac- 
celeration in the inner jets of blazars has been widely 
considered as a possible solution to this problem 
([Mannheim & Bicrmann 1992; Bcrczins ky et al.l 120061 : 
iDermer et al.1120091 : IMurase fc Ta kami 2009). In this pa- 
per, we argue that the detection of VHE radiation from 
the FSRQ 4C -1-21.35 supports this scenario, and pro- 
vides a solution to the rapid VHE variability and large 
7-ray luminosity as an effect of ultra- high energy proton 
accelerated in the inner jetI3 

A new feature of this study is the derivation of beam- 
ing factors for secondaries formed by photohadronic pro- 
cesses. Synchrotron and SSC emissions, and radiations 
from photohadronic processes with internal isotropic ra- 
diation flelds, have a beaming factor oc J^ due to com- 
pression in solid angle and time and enhancement in en- 

^ FSRQ blazars are scarce within the GZK volume, and their 
average emissivity is much less than blazar BL Lac objects, which 
may therefore be the preferred blaz ar class to accelerate most 
of the highest en ergy cosmic rays UDermer fc Hazzaguel 120101 : 
IMurase et~alll2012l) . 



ergy. External Compton scattering has a 5^ beaming fac- 
tor because Compton scattering is proportional to target 
photon energy density, which is boosted by two powers 
of fo when transformed to the comoving frame. Photo- 
hadronic processes with external isotropic photon sources 
have a beaming factor ex Jp, because photohadronic pro- 
duction is proportional to the target photon number den- 
sity, which is increased by a single power of the Dopper 
factor in the comoving frame. Threshold and spectral 
effects complicate the beaming factors, and the more de- 
tailed relations, including a formalism to calculate pro- 
duction spectra of photohadronic secondaries, are pre- 
sented in this paper. 

In Appendix 1^1 we also derived accurate photon pow- 
ers for bolometric 7-ray fluxes associated with syn- 
chrotron and SSC processes, and for external Comp- 
ton processes. These relations demonstrate the much 
stronger decline of flux for 7 rays from external Comp- 
ton scattering compared to SSC fluxes, which may help 
explain the relati ve number of misali gned radio galaxies 
of different types ()Abdo et al.ll2010c[ ). 

The main goal of this paper is to propose a new mecha- 
nism to generate rapidly variable 7 rays at large distances 
from a black-hole engine in order to explain the puzzling 
observations of 4C -1-21.35. The variability timescale 
is determined by processes in the inner jet that relate to 
the dynamical timescale of the black h ole. For blac k hole 
mass es of 1.5 - 8 x lO^Mg (|Wang et a l. 2004; Shaw et al.l 
|2012[) . the dynamical timescale « 1 - lOx larger than the 
variability timescale. Some of the same solutions applied 
to rapid variability in TeV bl azars, ie.z. iBegelman et aLl 
120081: iNaravan fc Piranll2012[ ). could operate in the inner 
jet of 4C -1-21.35 if the larger black-hole mass is correct. 

We assume that the inner jets of blazars accelerate 
UHECRs that undergo photohadronic losses with am- 
bient radiation field s and generate UHE neu trons, neu- 
trinos, and 7 rays (jAtovan fc Dermeil 120031 ). Outflow- 
ing UHECR neutrons undergo photohadronic losses with 
IR torus photons to make pions that decay into hyper- 
relativistic leptons. If the magnetic field is ~ 0.01 — 1 mC, 
as shown by the cross-hatched region in Figure ^ then 
the VHE synchrotron radiation preserves the rapid vari- 
ability of the inner engine at the pc scale. This mech- 
anism is proposed as the origin of the rapid variability 
of 70 - 400 GeV 7 rays observed with the MAGIC tele- 
scopes. Only a very narrow angular range of electrons 
will contribute to the observed VHE radiation, so that 
the variability of the central source is preserved even 
if the opening angle of the relativistic jet is not nar- 
row. UHE photons escaping to larger distances, e.g., 
at the Mpc-scale region in the large scale structure sur- 
rounding the source may furthermore produce a slowly 
variable VHE s ynchrotron pa ir echo with timescales of 
- 0.1-1 yr (|Murasel[2012l ). while UHE protons es- 
caping to the magnetized regions can give almost non- 
varia ble signals (jGabici fc Aharonianll2005t iKotera et al.l 
|2011|) . These radiations are useful as characteristic sig- 
nals of UHECR acceleration, and are potentially observ- 
able with the planned CTA (Ac tis et al. 2011). 

The 7-ray fluxes made principally by the synchrotron 
mechanism become less variable with decreasing energy, 
but at sufficiently low energies, 7 rays from the inner 
jet can avoid 77 absorption and additionally contribute 
to the observed flux. Indeed, Figure |8] shows that the 
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77 absorption for inner jet 7 rays becomes significant 
only above a few GeV. Tliis corresponds precisely to 
the range of energies where the Fermi-LAT has dis- 
covered spectral breaks in FSRQs a nd low- and inter- 
mediate synchrotron peaked blazars (jAbdo et al.l 120091 
I2010airbl ). Whether this is a coincidence, or a cause of 
the GeV spectral breaks, it is pertinent that most mod- 
els for the GeV br eaks invoke pro duction with i n the 
BLR (iPoutanen fc Stern 2010: Fi nke fc Dermeii 120101: 
lAckermann et al.ll2010t iStern fc Poutanenll2011| ). An in- 
ner jet origin of the GeV radiation is consistent with 
the model proposed here, and the synchrotron 7 rays 
at multi-GeV/TeV energies would appear as a separate 
spectral component emerging from the attenuated in- 
ner jet radiation spectrum, which furthermore display 
increased variability with energy. 

An accurate fit to the combined Fermi and MAGIC 7- 
ray spectrum requires detailed modeling of the radiation 
in the inner jet in addition to the subsequent cascade 
radiations formed as the VHE synchrotron photons pass 
through the IR photon field. Besides the leptonic emis- 
sions in the inner jet, a model involving the UHECR pro- 
tons that are accelerated and undergo photopion losses 
to make escaping UHECR neutrons must be included. 
As noted earlier, long-term average apparent 7-ray lumi- 
nosities > 10^® erg s~^ are common in FSRQs, so ap- 
parent jet powers must even be larger. The apparent 
luminosity in UHECRs accelerated in the inner jet of 4C 
+21.35 must be sa 10^^ - 10^^ erg s~^ to compensate for 
the inefficiencies to produce neutrons in the inner jet and 
pions at the pc scale. For the target photons from the 
external BLR and IR torus fields, the radiative efficiency 
for p7 photopion processes and the requirement of rapid 
variability from a single blob leads to the requirement 
of large Doppler factors, ^ 100 in order to make an ap- 
parent luminosity « 10^^ erg s~^ in outflowing neutrons 
from a jet whose absolute power is Eddington-limited. 

Internal shocks formed by colliding plasma shells 
ejected from the nuclear black hole are often considered 
as the mechanism that dissipates energy into broadband 
nonthermal radiation. The constraint R < r(5Dcivar/(l + 
z) on variable emission from colliding shells means that 
R < 0.01 fr/50)^(tv .. 7600 s) PC for F « <5d (e.g., 
iTavecchio et al.l 120101 ) . If the location of the emission 
site is required to be at a distance larger than a few pc 
from the central black hole, then unless F > 10"^, which 
exceeds by an order of magnitude the outflow Lorentz 
factor inferred from any blazar data, such a model for 
the MAGIC observations cannot explain the VHE radi- 
ation due to the severe attenuation. 

The necessity to reconcile the contrary features of 
short term variability, which points to an inner jet ori- 
gin, and detection of VHE 7 rays that must originate 
at the pc scale , has l ed to several proposed solutions. 
ITavecchio et al.l ()2011| ) consider a system with a com- 
pact zone outside the BLR and a second, possibly more 
extended zone either inside or outside the BLR to repro- 
duce the complete SED. The compact emitting regions 
making the VH E radiation must have large Doppler fac- 
tors, (5d ~ 75. iNalewaiko et al.l ()2012D consider various 
energetic constraints on 4C -1-21.35 system, and also ar- 
rive at the requirement of compact emitting zones at pc 
scales. The question then is the origin of these compact 



regions. One poss ibility is the formation of compact rec- 
ollim ation shocks ()MarscheHl2006l : iBromberg fc LevinsonI 
120091) that are produced when the external medium pres- 
sure overcomes the radiative jet, as proposed for M87 
flares or the radio cores of blazars. The large powers 
from 4C -1-21.35 challenge such an origin of the strongly 
variable 7-ray fluxes. These behaviors might also be 
reconciled in highly magnetized Foynting jet models, 
mini- j et models, or turbul ent ce lls (e.g., [G iannio s et al.l 
120091: iMarscher fc Jorstadl l2010t JNalewaiko et al.ll2011| l 
that have been considered in TeV BL Lac objects. 
INalewaiko et al.l (|2012| ) conclude that self-collimating jet 
structures might produce the conditions needed to explan 
the observations. 

Here we propose a new technique to produce rapid vari- 
ability of VHE 7 rays far from the central black hole 
as a result of UHECR neutron production in the inner 
jet, with the outflowing neutrons making cascade syn- 
chrotron radiation from leptons formed as secondaries 
from photopion interactions of the UHECR neutrons 
with the dust IR radiation field at the pc scale. The 
power in escaping neutrons is accompanied by 7-ray and 
neutrino radiations from the decay secondaries of neu- 
tral and charged pions. This approach to variability, 
which depends on rectilinear propagation of the parti- 
cles to large dis t ances , has similarities with an idea of 
iGhisellini et al.l ()2009D where rapid variability of TeV 
blazars like PKS 2155-304 or Mkn 501 results from elec- 
trons forced to follow field lines. This magnetocentrifugal 
model would not likely work with 4C -1-21.35 at the pc 
scale, but also involves a geometry essentially different 
from a relativistic plasma. 

This model predicts neutrino fluxes detectable with 
IceCube from flares of FSRQs, not only PKS 1222-f216, 
but also sources like 3C 279 and PKS 1510-089 that gen- 
erate a VHE fluence reaching « lO^'^ erg cm~^. For the 
one-half hour of activity of 4C +21.35, the total electro- 
magnetic fluence is w 10^°^ erg cm"'^, but the flare itself 
could have lasted for a much longer time than the half 
hour during which MAGIC was observing. VIoreover, the 
long-term jet radiation could accelerate UHECRs in the 
inner jet with associated neutrino production from FS- 
RQs with VHE emission. Furthermore, UHECR protons 
and ions could escape upstream to produce emissions at 
the pc scale without significant 7-ray and neutrino pro- 
duction in the inner jet, but would be strongly depleted 
by interactions in the dust torus. 

In conclusion, the proposed model, if correct, would 
provide an important clue to the question of the origin of 
UHECRs. The feature of transport of inner jet energy to 
the pc scale via UHECRs would explain why the question 
of the location of the 7-ray emission site in blazars has 
been so puzzling. UHECR processes could be essential in 
fitting the SED of all blazars, not just the FSRQs from 
which VHE radiation has been detected. 
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APPENDIX 
DERIVATION OF THE PHOTON POWER 
The apparent isotropic bolometric synchrotron luminosity is related to L'-^^ the isotropic synchrotron luminosity in 
the comoving jet frame, according to the relation Lsyn(^) = 5^L'^y^{Q,') ~ (5o(igy,-,/47r) = S^L'^y^/Ai:. The absolute 
synchrotron power radiated in all directions is, after multiplying by a factor of 2 for a two-sided jet, given by 

(Al) 

J \ 4:71 J J_l 

This can also be written as 

^syn.abs — j ^syn — „ c4 ^syn ' „ c4 ^syn j \-^^ ) 

noting Lgyn = 4:TrLsyn{^)- This beaming factor also applies to SSC emissions in the standard one- zone blazar model. 
The apparent isotropic bolometric Compton luminosity resulting from an isotropic comoving electron distribution 
scattering an external isotropic radiation field transforms according to Leg oa SqL'^q. Following the same reasoning, 
therefore, the absolute Compton power radiated in all directions is given in terms of the apparent Compton luminosity 
by the relation 

I\d^iS|, 2r4(5 + io/?^ + /?4) ^ /^^i32r4 

i-ECabs - ^e ^EC - ^^g -f-EC ^ "s^^^C • (Adj 

These relations apply to bolometric luminosities. Measurements over a finite frequency range introduce integration 
limits that give corrections that depend on spectral parameters. 

UHE NEUTRON BEAM PRODUCTION IN INNER JETS 

UHECR acceleration may occur in inner jets via shock acceleration or possibly magnetic reconnection. For tvar — 
600 s, the comoving size of the emission region is r'f, — (5Dcivar/(l + 2^) — 3.6 x 10 ((5D/20)(ivar/600 s)/(l -f z) cm. As 
long as we consider blob radii smaller than the radius of the BLR region, the photon field in inner jets is important. 
Though the nonthermal photon field in inner jet is uncertain, following Tavecchio et al. (2010), we take the maximum 
allowed photon field as L^ = 10^^'^ erg s~^ at Vf, = lO^'' Hz and a photon index (3 — 2.2 at i^ > i^;, for demonstration. 
The effective photon energy density at Vb in the inner jets is 

u'' 
rieff = -^ ^ 1.2 X 10^3 cm-3i^^y^^_4j. 5(^0/20) 'V.^^Vivar/eOO s) '(1 + zf (Bl) 

Then, the photomeson production efficiency in the inner jets is (e.g.. lAtovan fc DermeH[2003l ; iMurase fc Beacomll201Qf ) 



~0.2(<5D/20)-V-i4(ivar/600 s)-^ L%,^{E / E^^f-^ (B2) 

where E''^ ~ 1.5 x 10^° eV [6^, / 2Qf v^l^. This implies that r/p^, - 8 x 10"^ at Ep - 10^^ eV, which is not so efficient. 
Note that the corresponding photomeson production efficiency due to external photon fields (during the dissipation) 
would be low for smaller values of (5d. In fact, for interactions with the dust photon field, we have 

77p^ - 3 X 10-3(uo/e04(i?/10''' cm), (B3) 

where R « Fr'b is the distance to the blob. 
If neutron production mainly occurs in such a compact region, we have the neutron conversion efficiency as 

77p^^„ « 2.5r,p^ ~ 0.48(,5D/20)-V,;i4(^var/600 s)~^ L\^^{E / E^t^ (B4) 

High-energy neutrons are depleted via conversion to protons or further photomeson production loss. The neutron 
absorption efficiency rjn-y^p is expected to be comparable to rjp^^n, so the critical energy (where rjn-y^p = 1) is 
estimated to be 

e: ~ 2.3 X 10^0 eV {Su/20f'/'^;}i'{U.r/600 sf'{L%,,f'^' (B5) 

As a result, the resulting neutron spectrum can be written as 

ElN„iE„) « ^^^^^^hUn^E^NpiEp) oc i?^-+^ (B6) 
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where we have assumed that the maximum proton energy is lower than Ef^ for the last expression. The Hillas condition 
gives E™'^'^ w 0.8E™'^^ ^ 10^^'^ eV, below which we expect E^Nn oc E^-^ for s ~ 2. However, in our cases, the neutron 
beam production in inner jets typically requires quite large cosmic-ray luminosities. On the other hand, one may also 
expect proton escape especially at the maximum energies. But details depend on the uncertain escape mechanism. 
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